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Measurements of midrapidity charged particle multiplicity distributions, dNch/dg, and midrapid¬ 
ity transverse-energy distributions, dEr/dr], are presented for a variety of collision systems and 
energies. Included are distributions for Au-|-Au collisions at = 200, 130, 62.4, 39, 27, 

19.6, 14.5, and 7.7 GeV, Cu-fCu collisions at — 200 and 62.4 GeV, Cu-I-Au collisions at 

,/sffff = 200 GeV, U-l-U collisions at = 193 GeV, d-|-Au collisions at = 200 GeV, 

®He-l-Au collisions at = 200 GeV, and p-\-p collisions at = 200 GeV. Gentrality- 

dependent distributions at midrapidity are presented in terms of the number of nucleon partici¬ 
pants, Apart, and the number of constituent quark participants, Nqp. For all A-\-A collisions down 
to =7.7 GeV, it is observed that the midrapidity data are better described by scaling with 

Nqp than scaling with Apart. Also presented are estimates of the Bjorken energy density, £bj, and 
the ratio of dET/drj to dNch/drj, the latter of which is seen to be constant as a function of centrality 
for all systems. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

Systematic measurements of the centrality dependence 
of transverse energy production and charged particle 
multiplicity at midrapidity provide excellent character¬ 
ization of the nuclear geometry of the reaction and are 
sensitive to the dynamics of the colliding system. For ex¬ 
ample, measurements of dNch/dr] and dET/drj in Au-|-Au 
collisions at y/s/f/f = 200 GeV and 130 GeV as a func¬ 
tion of centrality expressed as the number of participant 
nucleons, Apart) exhibit a nonlinear increase with increas¬ 
ing Apart- This has been explained by a two-component 
model proportional to a linear combination of the num¬ 
ber of collisions, Acoii, and Apart In a previous 

study by the PHENIX collaboration, measurements of 
dEr/dij and dNc^/dij for Au-|-Au collisions at 200, 130, 
and 62.4 GeV are presented along with comparisons to 
the results of several models @|. The models that were 
examined included HIJING a final state parton sat¬ 
uration model called EKRT ra, an initial state parton 
saturation model called KLN [2| , and a multiphase trans¬ 
port model called AMPT Q. The comparisons showed 
that most models could reproduce some of the features 
of the data, but most failed in describing all of the data 
with the HIJING and AMPT models best describing the 
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overall trends, including the nonlinear increase of dEp/dr] 
and dNch/dij as a function of Apart- 

It was also proposed that dNch/dt] is linearly propor¬ 
tional to the number of constituent-quark participants 
without a significant contribution from a hard scatter¬ 
ing component Q. Recently, the PHENIX collaboration 
at Brookhaven National Laboratory’s Relativistic Heavy 
Ion Gollider (RHIG), presented dEp/dr] distributions at 
midrapidity for Au-|-Au collisions at y/sffff = 200, 130, 
and 62.4 GeV, d-|-Au collisions at = 200 GeV, and 

p-\-p collisions at — 200 GeV Q. The data are 

better described by a model based upon the number of 
constituent quark participants than by the wounded nu¬ 
cleon model 0. Here, this study will be extended to in¬ 
clude both dEp/drj and dNch/dr] measurements at midra¬ 
pidity in Au-|-Au collisions down to = 7.7 GeV. 

This study will also examine the centrality dependence 
of dEp/dij and dN^h/dr] for smaller systems including 
Gu-I-Au, Cu-kCu, d-l-Au, and ^He-I-Au. 

Recent lattice quantum chromodynamics (QGD) cal¬ 
culations indicate that the transition from quark to 
hadronic matter is a crossover transition at h igh tempera¬ 
ture and small baryochemical potential, ps [l0| . At high 
values of /is and low temperatures, model calculations 
indicate the presence of a first-order phase transition and 
the possibility of a critical end point in the QGD phase 
diagram [HI- Relativistic heavy ion collisions serve as 
excellent probes of the QCD phase diagram [I^. The re¬ 
gion of the QCD phase diagram sampled by the collisions 
can be controlled by changing the beam energy. Lowering 
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the beam energy corresponds to raising the value of /is. 
From 2010 to 2014, RHIC executed a beam energy scan 
program to explore the QCD phase diagram, look for ev¬ 
idence of the phase boundaries and search for evidence 
of the critical end point. Presented here are dEx/dri and 
dNch/drj measurements from the beam energy scan as a 
function of centrality expressed as the number of nucleon 
participants, IVpart, from Au-|-Au collisions at -s/s^ = 
200, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV. 

Over the past 15 years, PHENIX has collected a com¬ 
prehensive dataset covering a wide variety of colliding nu¬ 
clei and collision energies, including the Au-|-Au collision 
beam energy scan mentioned above. Presented here will 
be charged particle multiplicity and transverse energy 
measurements from the following systems: Au-|-Au colli¬ 
sions at = 200, 130, 62.4, 39, 27, 19.6, 14.5, and 

7.7 GeV; Cu-I-Cu collisions at = 200 and 62.4 GeV; 
Gu-I-Au collisions at = 200 GeV; IH-U collisions 

= 193 GeV; ^He-I-Au collisions at = 

200 GeV; d-l-Au collisions at = 200 GeV; and p+p 

collisions at = 200 GeV. The results will be dis¬ 

cussed in the context of scaling with the number of par¬ 
ticipant nucleons (iVpart) and the number of participant 
quarks {Ngp). 

PHENIX has previously published charged parti¬ 
cle multiplicity distributions from Au-|-Au collisions at 
= 200 GeV i, Au-|-Au collisions at = 

130 GeV i [g, and Au-|-Au collisions at = 

19.6 GeV Q. PHENIX has also previously published 
transverse energy distributions from Au-|-Au collisions 
at = 200 GeV Q, Au-I-Au collisions at 

= 130 GeV g, Au-fAu collisions at = 62.4 

GeV Au-I-Au collisions at = 19.6 GeV Q, and 

minimum-bias distributions for d-l-Au and p+p collisions 
at = 200 GeV 0. Here the previously published 

PHENIX results will be presented along with data from 
the many new collision systems in a consistent format to 
facilitate comparisons. 

Similar measurements have been published by the 
other RHIC experiments. Charged particle multiplic¬ 
ity distributions have been published by BRAHMS for 
Au-I-Au collisions at = 200 and 130 GeV T5|. 

STAR for Au-I-Au collisions at y/s^ =130 GeV Tq . 
and PHOBOS for Au-I-Au collisions at y/s^ = 200, 
130, 62.4, 56, and 19.6 GeV along with Cu-|-Cu colli¬ 
sions at y/Sfjfj = 200 and 62.4 GeV, d-l-Au collisions at 
y/s^ = 200 GeV, and p+p collisions at y/s^ = 410 and 
200 GeV g. Transverse energy distributions have been 
published by STAR for Au-I-Au collisions at y/s^ = 
200 GeV g] . Presented here are many collision systems 
and energies that have not been previously published by 
PHENIX or the other RHIC experiments, especially for 
the transverse energy measurements. The first complete 
results on charged particle multiplicity and transverse 
energy from the RHIC beam energy scan program con¬ 
ducted from 2010 to 2014 are also included. 

This paper is organized as follows. The PHENIX de¬ 
tector and the methods used for centrality determination 


in each dataset will be described in Section [TTl The analy¬ 
sis of the data to measure dEx/dr] and dNch/drj including 
a description of estimates of the systematic uncertainties 
is described in Section EH The centrality dependent re¬ 
sults at midrapidity from the Au-I-Au beam energy scan 
in terms of Vpart are presented in Section IIVI A descrip¬ 
tion of the centrality dependent results at midrapidity 
for Cu-I-Cu and Cu-|-Au collisions in terms of Vpart are 
found in Section|Vl A description of the centrality depen¬ 
dent results at midrapidity for U-l-U collisions in terms 
of Vpart are found in Section |Vll Section IVHI contains a 
description of the centrality dependent results at midra¬ 
pidity of ^He-|-Au and d-|-Au collisions in terms of Vpart- 
A review all of the centrality dependent results in terms 
of Nqp is presented in Section IVHII Section |T2 contains 
a summary of the results. Data tables for all data sets 
are tabulated in the Appendix. 


II. THE PHENIX DETECTOR 


The PHENIX detector comprises two central spec¬ 
trometer arms, two muon spectrometer arms, and a set 
of forward detectors. All of the detector components and 
their performance is described elsewhere [l9l| . The analy¬ 
sis of charged particle multiplicity utilizes detectors in the 
central arm spectrometer [20| including the drift cham¬ 
ber (DC) and pad chamber 1 (PCI) detectors. The drift 
chambers are cylindrically shaped and located radially 
from 2.0 to 2.4 meters. The DC covers the pseudora¬ 
pidity region \r]\ <0.35 and 90° in azimuth for each arm. 
The DC has a resolution better than 150 pm in better 
than 2 mm in the z direction, and a two track separation 
better than 1.5 mm. The PCI detector is a multiwire 
proportional chamber mounted on the outer radius of 
the drift chamber at 2.5 m from the beam axis. PCI cov¬ 
ers the full central arm acceptance. PCI measures min¬ 
imum ionizing particles with an efficiency greater than 
99.5% with a position resolution of 1.7 mm by 3 mm and 
a two track separation of 4 cm. Reconstructed tracks 
from the drift chamber with an associated hit from PCI 
are counted as charged particle tracks in the multiplicity 
measurement. 

The analysis of transverse energy utilizes five of 
the lead-scintillator (PbSc) electromagnetic-calorimeter 
(EMCal) sectors in the central arm spectrometers [El| . 
Each calorimeter sector covers a pseudorapidity range 
of \r]\ <0.38 and subtends 22.5° in azimuth for a total 
azimuthal coverage of 112.5°. The front face of each sec¬ 
tor is located 5.1 m from the beam axis. Each sector 
contains 2592 PbSc towers arranged in a 36 x 72 ar¬ 
ray. Each tower has a 5.535 x 5.535 cm surface area 
and a thickness of 0.85 nuclear interaction lengths or 
18 radiation lengths. The PbSc EMCal energy resolu¬ 
tion has been measured using test beam electrons to be 
^ 2.1%, with a measured response pro¬ 

portional to the incident electron energy to within ±2% 
over the range 0.3 < Eg < 40 GeV. 
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For all data sets, a minimum-bias tri gger is provided 
by a pair of beam-beam counters (BBC) [^. Each BBC 
comprises 64 individual Cerenkov counters. Each BBC 
covers 27r azimuthally and a pseudorapidity range of 3.0 
< |? 7 | < 3.9. For p+p, d-|-Au, and ^He-|-Au collisions, 
an event is required to have at least one counter fire in 
each BBC. For all other collisions, at least two counters 
must fire in each BBC. The event vertex is reconstructed 
with a resolution of 2.0 cm in p+p collisions and 0.5 mm 
in central Au-l-Au collisions using the timing information 
from the BBCs. All events are required to have an event 
vertex within 20 cm of the center of the detector. 

Centrality determination in the original = 

200 GeV and = 130 GeV Au-f Au PHENIX anal¬ 

ysis is based upon the total charge deposited in the 
BBCs and the total energy deposited in the zero-degree 
calorimeters (ZDC) [l^. The ZDCs are a pair of hadronic 
calorimeters that cover the pseudorapidity range I 77 I > 6 . 
For subsequent data sets taken after 2002, only the BBC 
information is used for the centrality determination, in¬ 
cluding the following data sets: Cu-|-Au at y/s^ = 
200 GeV, Cu-bCu at y/s^ = 200 GeV, U-bU at 
= 193 GeV, ^He-bAu at y/s^ = 200 GeV, and d-l-Au 

y/s^ = 200 GeV. As the collision energy decreases, 
the width of the pseudorapidity distribution of produced 
particles becomes more narrow (^ . As a result, for ener¬ 
gies below y/s^ = 130 GeV, the acceptance of the ZDC 
is reduced, therefore only the BBC information is used 
for Au-I-Au collisions at y/s]^ = 62.4 and 39 GeV, and 
for Cu-bCu collisions at y/s^ = 62.4 GeV. Below y/F^ 
= 39 GeV, the BBC acceptance becomes sensitive to the 
presence of beam fragments, which affects the linear re¬ 
sponse of the BBC to the centrality. To avoid this non¬ 
linear response, the reaction-plane detector (RXNP) 
is used for the centrality determination for Au-|-Au col¬ 
lisions at y/Sj^j^ =7.7 GeV, which was taken during the 

2010 running period. The RXNP comprises two sets of 
plastic scintillators positioned on either side of the colli¬ 
sion vertex. Each RXNP detector is arranged in 12 az¬ 
imuthal segments separated into an inner and outer ring. 
The RXNP has an azimuthal coverage of 27r. The pseu¬ 
dorapidity coverage is 1.5 < \r]\ <2.8andl.0< |? 7 | <1.5 for 
the inner and outer ring, respectively. A 2 cm thick lead 
converter is located directly in front of the RXNP scin¬ 
tillators with respect to the collision region, which allows 
the RXNP to also measure contributions from neutral 
particles through conversion electrons. The RXNP is de¬ 
signed to measure the reaction plane angle, but it can 
also function well as a centrality detector, because the 
magnitude of the total charge measured by the RXNP is 
dependent on the centrality of the collision. In order to 
minimize contamination from beam fragments, only the 
outer ring of the RXNP is used for centrality determina¬ 
tion for Au-bAu collisions at y/s^ = 7.7 GeV. For the 

2011 data taking period and later when the Au-|-Au data 
sets at y/s^ = 27, 19.6, and 14.5 GeV were collected, 
the RXNP was removed in order to install a silicon ver¬ 
tex detector, which was being commissioned during this 


time. So, for these two data sets, the multiplicity of hits 
in the PCI detector were used to determine the central¬ 
ity. A summary of the centrality detectors used for each 
dataset is included in Table HI 

III. DATA ANALYSIS 

Table U provides a summary of the data sets used in 
this analysis. For Au-bAu collisions at y/s^ = 62.4 GeV, 
the dEx/dri analysis uses data taken in 2004 Q and the 
dNch/drj analysis uses data taken in 2010. The number 
of events are those events that pass the minimum-bias 
trigger condition for the dataset and have an event vertex 
within 20 cm of the center of the detector. 


TABLE I. Summary of the data sets used in this analysis. 


(GeV) 

System 

YiB&r N events 

Centrality 

Trigger eff. 

200 

Au-I-Au 

2002 

270 k 

BBC-bZDC 

93 ± 3% 

200 

Au-bAu 

2004 

133 M 

BBC-bZDC 

93 ± 3% 

130 

Au-bAu 

2000 

160 k 

BBC-bZDC 

93 ± 3% 

62.4 

Au-bAu 

2004 

20 M 

BBC 

86 ± 3% 

62.4 

Au-bAu 

2010 

12 M 

BBC 

86 ± 3% 

39 

Au-bAu 

2010 

132 M 

BBC 

86 ± 3% 

27 

Au-bAu 

2011 

24.5 M 

PCI 

86 ± 3% 

19.6 

Au-bAu 

2011 

6.3 M 

PCI 

86 ± 3% 

14.5 

Au-bAu 

2014 

6.8 M 

PCI 

85 ± 3% 

7.7 

Au-bAu 

2010 

803 k 

RXNP 

75 ± 3% 

200 

Cu-bCu 

2005 

558 M 

BBC 

93 ± 3% 

62.4 

Cu-bCu 

2005 

175 M 

BBC 

88 ± 3% 

200 

Cu-bAu 

2012 

2.6 B 

BBC 

93 ± 3% 

193 

U-bU 

2012 

317 M 

BBC 

93 ± 3% 

200 

®He-bAu 2014 

1.6 B 

BBC 

88 ± 4% 

200 

d-bAu 

2008 

1.4 B 

BBC 

88 ± 4% 

200 

P+P 

2003 

14.6 M 

— 

54.8 ±5.391 


A. Transverse Energy Analysis 

The analysis procedure for dEx /dr] is described in de¬ 
tail in Ref. Q and summarized here. The absolute en¬ 
ergy scale for each EMGal sector is calibrated using the 
7 r° mass peak from pairs of reconstructed EMGal clus¬ 
ters for each dataset. The transverse energy for each 
event was computed using clusters in the EMCal with an 
energy greater than 30 MeV composed of adjacent tow¬ 
ers each with a deposited energy of more than 10 MeV. 
Faulty towers and all towers in a 3 x 3 tower area around 
any faulty tower are excluded from the analysis. The 
transverse energy Ex is a multiparticle variable defined 
as the sum 

Ex = ^E^ sin 9^ 






















7 


dEr {r]) /dr] = sin 6 [t]) dE{T])/dr] , (1) 

where 0^ is the polar angle, ry = — lntan(0/2) is the 
pseudorapidity, Ei is by convention taken as the kinetic 
energy for baryons, the kinetic energy + 2 toat for an- 
tibaryons, and the total energy for all other particles, 
where rriN is the nucleon mass. The sum is taken over all 
particles emitted into a fixed solid angle for each event. 
An example of the raw Etemc distributions as a function 
of centrality for Au+Au collisions at = 14.5 GeV 

are shown in Fig. [TJa). 

In order to obtain the total hadronic Ex within a ref¬ 
erence acceptance of Ary = 1.0, At/) = 27r from the mea¬ 
sured raw transverse energy, AtemCi the total correction 
can be decomposed into three main components. First 
is a correction by a factor of 4.188 to account for the 
fiducial acceptance in azimuth and pseudorapidity. Sec¬ 
ond, a correction factor is applied to account for disabled 
calorimeter towers not used in the analysis. Third is a 
factor, k, which is the ratio of the total hadronic Ex in 
the fiducial aperture to the measured Exemc ■ Details on 
the estimate of the values of the k factor are given below. 

The k factor comprises three components. The first 
component, denoted krespouse, is due to the fact that 
the EMCal was designed for the detection of electromag¬ 
netic particles M- Hadronic particles passing through 
the EMCal only deposit a fraction of their total energy. 
The average EMCal response is estimated for the various 
particle species using the HIJING Q event generator for 
y/s^ above 7.7 GeV and the URQMD [23 event gen¬ 
erator for Au-I-Au collisions at y/s^ =7.7 GeV. The 
event generator output is processed through a GEANT- 
based Monte Carlo simulation of the PHENIX detector. 
For all of the data sets, 75% of the total energy incident 
on the EMCal is measured, thus k^esponse = 1/0.75 = 
1.33. The second component of the k factor, denoted 
kin flow, is a correction for energy inflow from outside 
the fiducial aperture of the EMCal. This energy inflow 
has two sources: from parent particles with an original 
trajectory outside of the fiducial aperture whose decay 
products are incident within the fiducial aperture, and 
from particles that reflect off of the PHENIX magnet 
poles into the EMCal fiducial aperture. The energy in¬ 
flow contribution is 24% of the measured energy, thus 
kinfiow = 1-0.24 = 0.76. The third component of the 
k factor, denoted kiosses, is due to energy losses. There 
are three components to the energy loss: from particles 
with an original trajectory inside the fiducial aperture of 
the EMCal whose decay products are outside of the fidu¬ 
cial aperture (10%), from energy losses at the edges of 
the EMCal (6%), and from energy losses due to the en¬ 
ergy thresholds (6%). The total contribution from energy 
losses is 22%, thus kiosses = l/(l-0.22) = 1.282. The to¬ 
tal k factor correction is fc = Response X hnflow X kiosses 
= 1.30. This value varies by less than 1% for all data 
sets. 

There are several contributions to the systematic un¬ 
certainties for the dEx/dr] measurement which are added 
in quadrature to obtain the total uncertainty. These 


contributions include the following: uncertainties due to 
the energy response of the EMCal, uncertainties due to 
the estimate of the EMCal acceptance, uncertainties due 
to the estimate of losses and inflow, uncertainties due 
to sector-by-sector variations, uncertainties due to the 
noise background estimate, uncertainties due to the trig¬ 
ger background estimate, and uncertainties due to the 
trigger efficiency estimate. A summary of the systematic 
uncertainties for the dEx/dr] analysis of each dataset is 
listed in Table mi for each dataset and further explained 
below. 

There is an uncertainty due to the energy response 
of the EMCal. This includes uncertainties in the ab¬ 
solute energy scale, uncertainties in the estimate of the 
hadronic response, uncertainties from energy losses on 
the EMCal edges and uncertainties from energy thresh¬ 
olds. The uncertainties in the hadronic response include 
a 3% uncertainty estimated using a comparison of the 
simulated energy deposited by hadrons with different mo¬ 
menta with test beam data [2ll | along with an additional 
1% uncertainty in the particle composition and momen¬ 
tum distribution. There is an estimated uncertainty of 
2% for the calculation of the EMCal acceptance. There is 
an estimated uncertainty of 3% for the calculation of the 
fraction of the total energy incident on the EMCal fidu¬ 
cial area (losses and inflow). There is an uncertainty due 
to sector-by-sector variations in the energy measurement. 
There is an uncertainty due to the noise, or background, 
contribution which is estimated to be consistent with zero 
with uncertainties determined by measuring the average 
energy deposited per sector in events where all the par¬ 
ticles are screened by the central magnet pole tips by 
requiring an interaction z-vertex of -1-50 < z < -|-60 cm 
and —60 < z < —50 cm. There is a centrality-dependent 
uncertainty for background due to multiple interactions 
and trigger effects. 

There is also an uncertainty in the trigger efficiency de¬ 
termination. The method by which the trigger efhciency 
is calculated is described in Q. The BBC trigger effi¬ 
ciency for Au-I-Au collisions ranges from 93% at 
= 200 GeV to 75% at = 7.7 GeV. The trigger 

efficiencies for each dataset are summarized in Table ID 
Note that the trigger inefficiency leads to a partial loss 
of the more peripheral collisions while the trigger is fully 
efficient for midcentral and central collisions. Because 
the centrality is defined for a given event as a percent¬ 
age of the total geometrical cross section, an uncertainty 
in the trigger efficiency translates into an uncertainty in 
the centrality definition. This uncertainty is estimated 
by measuring the variation in dEx/dr] by redefining the 
centrality using trigger efficiencies that vary by ±1 stan¬ 
dard deviation. 

The trigger efficiency uncertainty allows for bend¬ 
ing or inclination of the points. So, when plotting 
{dEx/dr])/{0.5Npart) and {dNeh/dr])/{0.5Npe,rt), the trig¬ 
ger efficiency will be represented by error bands about the 
points within which the points can be tilted. The other 
systematic and statistical uncertainties are represented 







= 14.5 GeV are shown in 




FIG. 1. (Color online) Raw Etemc (a) and Nch (b) distribu¬ 
tions for y^SNN = 14.5 GeV Au-I-Au collisions. Shown are the 
minimum-bias distribution along with the distributions in 5% 
wide centrality bins. All the plots are normalized so that the 
integral of the minimum-bias distribution is unity. 


by error bars. 


B. Charged Particle Multiplicity Analysis 

In previous PHENIX publications [1, [l^ for Au-|-Au 
collisions at = 200 and 130 GeV, charged parti¬ 

cle multiplicity was measured using cluster pairs recon¬ 
structed from the PCI and PCS detectors in the absence 
of a magnetic field. The dNch/dr] values quoted here for 
Au-I-Au collisions at = 200 and 130 GeV are from 

the previous analyses. For all other collision species and 
collision energies, charged particle multiplicity is mea¬ 
sured using reconstructed tracks from the drift chamber 
that have an unambiguous match to a reconstructed clus¬ 
ter in the PCI detector with the magnetic field turned 
on. In order to remove multiple counting of incorrectly 
reconstructed tracks in the drift chamber, commonly re¬ 
ferred to as ghost tracks, a charge-dependent track prox¬ 
imity cut is applied. The two methods give consistent 
results for 200 GeV Au-I-Au collisions. An example of 
the raw Nch distributions as a function of centrality for 


the Au-f An collisions at 
Fig. [11(b). 

In order to obtain the total charged particle A^ch within 
a reference acceptance of Arj = 1.0, A(^ = 27r from the 
measured raw multiplicity, five corrections are applied. 
First is a correction of 3.74 to account for the fiducial 
acceptance in azimuth and pseudorapidity. The second 
correction is applied to account for drift chamber and 
PCI inefficiencies within the fiducial acceptance. The 
third correction is applied to account for particles with 
a transverse momentum below the 200 MeV/c minimum 
Pt cut applied to reconstructed tracks. This correction 
is determined using the average of results from the HI- 
JING event generator Q and the URQMD event gen¬ 
erator [ 2 ^ to estimate the fraction of the total charged 
particle multiplicity lying below px = 200 MeV/c. The 
collision energy cutoff for the HIJING event generator 
lies above = 7.7 GeV, so only URQMD is used for 

Au-fAu collisions at ^/s^ = 7.7 GeV. This correction 
is 22% for Au-I-Au collisions at y/s^ = 62.4 GeV and 
23% for Au-bAu collisions at =7.7 GeV. There 

is an estimated 2% uncertainty for this correction. The 
fourth correction is a centrality-dependent correction for 
the track reconstruction efficiency. 

The last correction is an in-flight decay correction that 
accounts for particle decays after the collision interaction 
that can add or remove charged particles from the mea¬ 
sured multiplicity. This includes primary charged parti¬ 
cles that decay and miss the detector. It also includes 
feed-down from neutral primary particle decays that go 
into the detector. This correction is determined by pro¬ 
cessing simulated events from the HIJING Q event gen¬ 
erator for y/s^ above 7.7 GeV, and the URQMD 
event generator at y/s^ = 7.7 GeV. Below y/s^ = 
62.4 GeV, results from the two event generators are con¬ 
sistent with each other within the uncertainties. The 
event generator output is processed through a GEANT- 
based simulation of the PHENIX detector response. For 
Au-fAu collisions, this correction varies from 0.99 at 
y/s^ = 200 GeV to 1.061 at y/s^ = 7.7 GeV. The 
energy dependence is primarily due to the decrease of 
the particle momenta and the narrowing of the width of 
the T] distribution at lower energies that affects the num¬ 
ber of tracks from the decay of particles coming from 
comparable rapidities. 

There are several contributions to the systematic un¬ 
certainties for the dNch/drj measurement which are added 
in quadrature to obtain the total uncertainty. A sum¬ 
mary of the systematic uncertainties for the dNch/dr] 
analysis for all data sets is listed in Table IHII There is 
an estimated uncertainty of 4% for the acceptance cor¬ 
rection. There is an uncertainty for the estimate of the 
correction for in-flight decays that varies from 2.9% at 
y/s^ = 200 GeV to 5.9% at y/s^ = 7.7 GeV. There 
is a 2% uncertainty for the estimate of charged particle 
multiplicity for low px below 200 MeV/c. There is a 
centrality dependent uncertainty due to the occupancy 
of the PGl detector that varies from 3.5% to 1.2% for 
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TABLE II. Summary of the systematic uncertainties for the dEx/di] measurement for each dataset, given in percent (%). If a 
range is specified, the value for central collisions is listed first and the value for the most peripheral collisions presented for the 
dataset is listed second. If no value is specified, then there is no contribution to the systematic uncertainty for that dataset. 


Dataset 

energy resp. 

acceptance 

losses & inflow 

sector-by-sector 

noise 

trigger bkg. 

trigger eff. 

200 GeV Au+Au 

3.9 

2.0 

3.0 

— 

0.2-6.0 

— 

0.3-16.0 

130 GeV Au+Au 

3.8 

2.0 

3.0 

— 

0.4-10.0 

— 

0.3-16.0 

62.4 GeV Au+Au 

4.3 

2.0 

3.0 

2.2 

0.4-4.1 

0.01-0.06 

0.3-16.1 

39 GeV Au+Au 

4.5 

2.0 

3.0 

1.6 

0.5-3.6 

0.002-0.02 

0.2-16.3 

27 GeV Au+Au 

4.5 

2.0 

3.0 

2.2 

0.5-3.5 

0.006-0.04 

0.3-13.1 

19.6 GeV Au+Au 

4.7 

2.0 

3.0 

2.8 

0.5-3.5 

0.008-0.07 

0.3-13.4 

14.5 GeV Au+Au 

4.7 

2.0 

3.0 

2.9 

0.5-3.4 

0.007-0.04 

0.3-9.8 

7.7 GeV Au+Au 

4.7 

2.0 

3.0 

3.7 

0.5-3.4 

0.002-0.05 

0.4-10.6 

200 GeV Cu+Cu 

3.9 

2.0 

3.0 

5.9 

0.2-6.0 

0.002-0.04 

0.3-6.5 

62.4 GeV Cu+Gu 

4.5 

2.0 

3.0 

2.2 

0.4-4.1 

0.006-0.02 

0.3-8.1 

200 GeV Cu+Au 

3.9 

2.0 

3.0 

2.8 

0.5-3.5 

0.02-0.20 

0.2-8.8 

193 GeV U+U 

3.9 

2.0 

3.0 

2.5 

0.2-6.0 

0.001-0.03 

0.4-9.3 

200 GeV d+Au 

3.9 

2.0 

3.0 

6.5 

0.2-0.2 

0.13-0.21 

0.3-5.1 

200 GeV ®He+Au 

3.9 

2.0 

3.0 

3.9 

0.2-0.2 

0.08-0.16 

0.2-5.2 

200 GeV p+p 

3.9 

2.0 

3.0 

3.9 

0.2 

0.60 

— 


TABLE III. Summary of the systematic uncertainties for charged particle multiplicity for each dataset given in percent (%). If 
a range is specified, the value for central collisions is listed first and the value for the most peripheral collisions presented for 
the dataset is listed second. If no value is specified, then there is no contribution to the systematic uncertainty for that dataset. 


Dataset 

acceptance 

decays 

low pt 

occupancy 

tracking eff. 

trigger bkg. 

trigger eff. 

200 GeV Au+Au 

2.3 

2.9 

2.0 

3.5-0.10 

— 

1.0 

0.3-16.0 

130 GeV Au+Au 

2.5 

2.5 

2.0 

3.1-0.10 

— 

1.0 

0.3-16.0 

62.4 GeV Au+Au 

4.0 

5.0 

2.0 

3.5-0.10 

5.0 

0.001-0.03 

0.2-16.1 

39 GeV Au+Au 

4.0 

5.4 

2.0 

3.0-0.03 

5.0 

0.001-0.009 

0.2-13.0 

27 GeV Au+Au 

4.0 

5.6 

2.0 

2.0-0.01 

5.0 

0.01-0.03 

0.2-13.3 

19.6 GeV Au+Au 

4.0 

5.7 

2.0 

1.9-0.01 

5.0 

0.002-0.003 

0.2-9.3 

14.5 GeV Au+Au 

4.0 

5.8 

2.0 

1.9-0.01 

5.0 

0.001-0.007 

0.3-9.8 

7.7 GeV Au+Au 

4.0 

5.9 

2.0 

1.2-0.01 

5.0 

0.001-0.03 

0.4-12.3 

200 GeV Cu+Cu 

4.0 

2.9 

2.0 

1.5-0.01 

5.0 

0.03-0.08 

0.3-8.0 

62.4 GeV Cu+Cu 

4.0 

5.0 

2.0 

1.0-0.01 

5.0 

0.02-0.01 

0.3-9.2 

200 GeV Cu+Au 

4.0 

2.9 

2.0 

2.6-0.05 

5.0 

0.001-0.07 

0.9-10.1 

193 GeV U+U 

4.0 

2.9 

2.0 

3.5-0.10 

5.0 

0.001-0.01 

0.4-9.3 

200 GeV d+Au 

4.0 

2.9 

2.0 

0.1-0.01 

5.0 

0.001-0.001 

0.3-7.2 

200 GeV ^He+Au 

4.0 

2.9 

2.0 

0.1-0.01 

5.0 

0.001-0.001 

0.2-6.5 

200 GeV p+p 

4.0 

2.9 

2.0 

0.01 

5.0 

0.0015 

— 


Au+Au central collisions from = 200 to 7.7 GeV. 

There is an estimated 5% uncertainty for the tracking ef¬ 
ficiency estimate. There is a centrality-dependent uncer¬ 
tainty for background due to trigger effects and multiple 
interactions. Finally, there is an uncertainty for the de¬ 
termination of the trigger efficiency, which is estimated 
in the same manner as for the dEx/dri analysis. 


IV. AU-bAU BEAM ENERGY SCAN RESULTS 

This section will present dEx/dr] and dNch/drj mea¬ 
surements as a function of centrality expressed as the 
number of nucleon participants, Ypart, from the RHIC 
beam energy scan that includes Au-|-Au collisions at 
= 200, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV. A 
Monte Carlo Glauber model calculation is used to obtain 
estimates of A^part as a function of centrality using the 
procedure outlined in Ref. [1^. At each collision energy, 
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FIG. 2. (Color online) {(IEt/ drj)/{Q.^Nps.T\.) (a) and {dNch/dr})/{0.5Npi^rt) (b) at midrapidity as a function of A^part for Au+Au 
collisions 200, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV. The lines bounding the points represent the trigger efficiency uncertainty 
within which the points can be tilted. The error bars represent the remaining total statistical and systematic uncertainty. 
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FIG. 3. The ratio of = 200 GeV Au+Au collisions to — 7.7 GeV Au+Au collisions for {dEx/drj)/{0.5Nps,rt) (a) 

and for (dAch/d»7)/(0.5Apart) (b). The lines bounding the points represent the trigger efficiency uncertainty within which the 
points can be tilted. The error bars represent the total statistical and systematic uncertainties. 


the Glauber model is run using the inelastic nucleon- 
nucleon cross sections, cr™^^ listed in Table HVl 

When plotting dET/dij and dNch/drj, systematic un¬ 
certainties are decomposed into two types. Type A un¬ 
certainties include point-to-point uncertainties that are 
uncorrelated between bins and include only statistical un¬ 
certainties in this analysis. The remaining uncertainties 
are classified as Type B uncertainties that are correlated 
bin-by-bin such that all points move in the same direc¬ 
tion, but not necessarily by the same factor. Because the 
magnitudes of the Type A statistical uncertainties are 
small compared to the magnitudes of the Type B uncer¬ 
tainties, the error bars in the plots presented below will 


represent the total statistical and systematic uncertain¬ 
ties added in quadrature. The trigger efficiency uncer¬ 
tainty is represented separately by error bands bounding 
the points within which the points can be tilted, as de¬ 
scribed in Section uni 

Examining the A^part dependence of dEx/dr] and 
dNch/drj normalized by the number of nucleon partici¬ 
pant pairs at midrapidity is useful to determine if the 
data scales by A^part and if the scaling changes as a func¬ 
tion of ^/s^. The results for Au+Au collisions for all 
beam energies at midrapidity are shown in Fig. [2] as a 
function of Afpart- For all energies, dET/dr] and dNc^i/di] 
do not scale with A^part; the magnitudes of dEx/dr] and 
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FIG. 4. (Color online) The excitation function of (dEr/dri)/{0.5 Np^r-t) (a) and (dA^ch/d77)/(0.5A^part) (b) for central col¬ 
lisions at midrapidity as a function of The error bars represent the total statistical and systematic uncertainties. 

For {dEr / duf) / {0.5NpB.Tt) (a), data are shown from FOPI [SJ, E80 2 [^. NA49 [H, STA R and CMS [^l^. For 
(dNch/dri)/(0.5Npg.rt) (b), data are shown from FOPI [2^, E802 (^. IsiLl^. NA49 [^, STAR [l§! |33|| . PHOBOS [i3i AL¬ 
ICE 0, and ATLAS ll|. 


TABLE IV. 

yJ^NN ■ 

Summary of the 

cross sections as 

a function of 

Energy 

[mb] 

[mb] 

47' [mb] 

200 

52.5 

42.3 

8.17 

130 

48.7 

39.6 

7.54 

62.4 

43.6 

36.0 

6.56 

39 

41.2 

34.3 

6.15 

27 

39.8 

33.2 

5.86 

19.6 

39.0 

32.5 

5.70 

15.0 

38.5 

32.0 

5.58 

7.7 

38.6 

31.2 

5.35 


dNch/drj increase as iVpart increases. It has been previ¬ 
ously observed that the shape of the distributions as a 
function of iVpart are preserved in Au-|-Au collisions from 
= 200 GeV to = 19.6 GeV dill. FigureUi 
shows the ratio of {dET/dr])/{0.5Npart) from Au-|-Au col¬ 
lisions at = 200 GeV to = 7.7 GeV, illus¬ 

trating that the shapes of the distributions are preserved 
down to = 7.7 GeV. FigurelSh shows the same for 

((iVch/d??)/(0.5Vpart). Previous measurements in fixed 
target h-|-A collisions showed that the total charged par¬ 
ticle multiplicity does scale well as a function of Vpart 
in the range of 10 <^/s^ < 20 GeV [36|. However, 
this measurement was made over the full rapidity range 
rather than at midrapidity. For the midrapidity measure¬ 
ments presented here, the Vpart scaling behavior does 


not change significantly from = 200 GeV down to 

= 7.7 GeV. 

Excitation functions of {dEr / dr]) / {0.5 Np^^rt) and 
(dVch/d??)/(0.5Vpart) are shown in Fig. S) Shown are 
the PHENIX data along with results from other experi¬ 
ments. The data points for the lower energies are from 
estimates described in Ref. i. Eor (dE;T/d7?)/(0.5Vpart), 
data are shown from FOPI 0-1% centrality Au-I-Au 
collisions [l^, E802 0%-5% centrality Au+Au colli¬ 
sions [ 2 ^, NA49 0-7% centrality Pb+Pb collisions 1^ . 

STAR 0%-5% centrality Au-fAu collisions [l8j| . 
and CMS 0%-5% centrality Pb-t-Pb collisions [s^. For 
(dNrh /dn) /( 0. SfVpart, ), data are shown from FOPI [T7I |. 
E802 Hi, [Hi, m, NA49 (H, STArH, ill, PHO¬ 
BOS 0-3% centrality Au-|-Au collisions [l7j| . ALICE 0%- 
5% centrality Pb-|-Pb collisions [s^, and ATLAS [s^ 
Pb-|-Pb collisions interpolated to 0%-5% centrality. The 
data are plotted on a log-log scale to illustrate the 
power law behavior of both {dET/dri)/{0.5Npait) and 
{dNch/dr])/{0.5Np^rt) as a function of \og{,/s):)^) for 
at or above 7.7 GeV. For {dET/dr])/{0.5Nps,rt), 
the data between yfs^ = 7.7 and 200 GeV are de¬ 
scribed by (dF;T/d7?)/(0.5Vpart)(ys)^) cx ^ 

where b = 0.428 ± 0.021. For {dNch/dr])/{0.5Npart), the 
data between = 7.7 and 200 GeV are described 

by {dNch/dr])/{0.5Npa.rt)iy/s):^) oc where 

b = 0.374 ± 0.028. The data deviate from the power 
law behavior below the lowest PHENIX measurement at 
^ = 7.7 GeV. 

The ratio of dEx/dr] to dNch/dr], referred to here sim- 
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FIG. 5. The Er/Nch ratio as a function of A^part for Au+Au collisions at varying values of 44 ^. The lines bounding the 
points represent the trigger efficiency uncertainty within which the points can be tilted. The error bars represent the total 
statistical and systematic uncertainties. 
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FIG. 6. (Color online) The Ex/Nch ratio as a function of 
for central Au+Au collisions and Pb+Pb collisions at 
midrapidity. The error bars represent the total statistical and 
systematic uncertainties. The Large-Hadron-Collider, LHC, 
data point has been obtained by taking the ratio of the CMS 
dEr/dr] data with the average of the ALICE and 
ATLAS [ 35 II data. For (dEr/dr))/(0.5Nua,Tt), data are taken 
from FOPI [13], E802 NA49 [H, [3|, STAR [3, and 
CMS For (dNch/ dv)J {0-5N^a.Tt), data are taken from 

FOPI 113, E802 E1J3,[33|, NA49 & STAR [3.113. PHO- 
BOS [i3; ALICE[^and ATLAS [sl. 


ply as Ex/Nch, is a variable that is related to the av¬ 
erage transverse mass of the produced particles [!]■ In 
previous measurements, this ratio has been observed to 
be independent of centrality and independent of 
in Au-I-Au collisions from = 200 to 19.6 GeV Q- 

Figure [5] plots the Ex/Nch ratio as a function of iVpart 
for Au-I-Au collisions at various values of ^/s^. For all 
cases, the ratio is constant with A^p^rt within the system¬ 
atic uncertainties. The excitation function of Ex/Nch 
is shown in Fig. [6] Here, the Large-Hadron-Collider 
point has been obtained by taking the ratio of the CMS 
dEx/d-q dat a with the average of the ALICE [s^ 
and ATLAS [S^l data. The ratio increases below y/s^ 
~ 10 GeV, levels off, and then increases at = 

200 GeV. 

The energy density per unit volume in nuclear colli¬ 
sions can be estimated from the energy density per unit 
rapidity (33. The Bjorken energy density can be calcu¬ 
lated as follows: 


£bj 


1 

A±t 




dEj' 

dr] 


( 2 ) 


where Aj_ is the transverse overlap area of the nuclei 
determined from the Glauber model, r is the formation 
time, and J(y, ry) is the Jacobian factor for converting 
pseudorapidity to rapidity. 


The Jacobian factor depends on the momentum distri¬ 
butions of the produced particles, which are dependent 
on the beam energy. The Jacobian factor for each beam 
energy in the PHENIX acceptance has been estimated us¬ 
ing the URQMD event generator, which well reproduces 
measured particle spectra over the RHIG beam energy 
range and, unlike HIJING, is valid at y/s^ = 7.7 GeV. 
Calculations of the Jacobian factor using URQMD are 
consistent with previous calculations using the HIJING 
event generator [3| ■ There is an estimated uncertainty of 
3% for this calculation for all beam energies. The values 
of the Jacobian factors are summarized in Table El 


TABLE V. Summary of the Jacobian scale factor estimated 
for each beam energy. 


Dataset 

J(y,6) 

200 GeV Au-I-Au 

1.25 

130 GeV Au-I-Au 

1.25 

62.4 GeV Au+Au 

1.25 

39 GeV Au+Au 

1.27 

27 GeV Au+Au 

1.27 

19.6 GeV Au+Au 

1.28 

14.5 GeV Au+Au 

1.30 

7.7 GeV Au+Au 

1.35 


The transverse overlap area is estimated using the a 
Monte Carlo Glauber model as A± ^ <Jx<^y, where Gx and 
Gy are the widths of the x and y position distributions 
of the participating nucleons in the transverse plane. A 
normalization to 'kB? , where R is the sum of the radius 
and a surface diffuseness parameters of the Woods-Saxon 
parametrization 


p(r) = l/(l + e(^-’-'‘)/“), (3) 

of the nuclear density profile, p(r), was applied for the 
most central collisions at impact parameter 6 = 0. 

A compilation of the Bjorken energy density multiplied 
by T for Au-|-Au collisions at various collision energies is 
shown in Fig. [T] The value of £bj increases with increas¬ 
ing yfs^ and also with increasing TVpart- The value of 
£bj for the most central Au-fAu collisions at y/s^ = 
7.7 GeV is 1.36 ± 0.14, which is still above the value of 
1.0 for a formation time of 1 fm/c that had been the pro¬ 
posed value above which the Quark-Gluon Plasma can be 
formed in Bjorken’s original paper [ 33 . It is also above 
the result of 0.7 ± 0.3 GeV/fm^ for the critical energy 
density obtained from lattice QCD calculations [3, 1^ - 
The excitation function of £bj multiplied by r is shown 
in Fig. 3 The results are shown on a log-log scale to illus¬ 
trate that £bj follows a power law behavior from y/s^ = 
7.7 GeV up to = 2760 GeV, ebjt oc 

where 6 = 0.422 ± 0.035. 
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FIG. 7. (Color online) The Bjorken energy density, ebj, mnl- 
tiplied by r as a function of A^part for An+Au collisions at 
varying values of The error bars represent the total 

statistical and systematic uncertainties. 


10 


o 

E 

a> 

O 


• PHENIX 0-5%Au+Au 
A CMS 0-5% Pb+Pb 


10 


10 ^ 


ils^N [GeV] 


10 " 


FIG. 8. (Color online) The Bjorken energy density, £bj, 
multiplied by r as a function of for central An+Au 

(PHENIX) and Pb+Pb (CMS) [3Ci collisions at midrapidity. 
The error bars represent the total statistical and systematic 


uncertainties. 




FIG. 9. (Color online) (di?T/dr?)/(0.5Apart) (a) and (dAch/d77)/(0.5Apart) (b) at midrapidity as a function of 
A^part for Cn+Cu and Cu+Au collisions. Also shown are results from Au+Au collisions at = 200 GeV for 

comparison. The lines bounding the points represent the trigger efficiency uncertainty within which the points can 
be tilted. The error bars represent the remaining total statistical and systematic uncertainty. 


V. RESULTS FOR CU+AU AND CU+CU 
COLLISIONS 

Measurements of dNch/drj in systems lighter than 
Au have been pnblished by PHOBOS for 200 GeV 
and 62.4 GeV Cu+Gu collisions El, showing that the 
Gu+Cu dNch/dr] distribution as a function of Apart ex¬ 


hibits similar features when compared to Au+Au colli¬ 
sions. Here, those measurements are extended to include 
measurements of dET/drj and the addition of measure¬ 
ments from the asymmetric Cu+Au system at = 

200 GeV. 

Figure [S] shows {dEx / drj) / and 

((iAch/di7)/(0.5Apart) at midrapidity as a function 
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FIG. 11. (Color online) The Bjorken energy density, ebj, 
mnltiplied by r as a fnnction of A^part for Cu+Cn, Cu+Au, 
and Au+Au collisions. The error bars represent the total 
statistical and systematic uncertainties. 


FIG. 10. The Er/Nch ratio as a function of A^part for Cu+Au 
collisions at = 200 GeV (a), Cu+Cu collisions at ^/Sjy^ 

= 200 GeV (b), and Cu+Cu collisions at = 62.4 GeV 

(c). The lines bounding the points represent the trigger ef¬ 
ficiency uncertainty within which the points can be tilted. 
The error bars represent the total statistical and systematic 
uncertainties. 


of A^part for Cu+Cu and Cu+Au collisions. Also 
shown for comparison are the data for Au+Au col¬ 
lisions at ^/s^ = 200 GeV. Both plots exhibit the 
trend established in Au+Au collisions of increasing 
{dET/dri)/{0.5Npi,rt) and {dNch/dr])/{0.5Nps.rt) with 


increasing Afpart and increasing The Cu+Cu 

and Cu+Au distributions at = 200 GeV are 

consistent with each other within the uncertainties of 
the measurement. All of the species (Au+Au, Cu+Au, 
and Cu+Cu) at y/s^ = 200 GeV are consistent with 
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each other for all overlapping values of A^part- This 
behavior had been previously noted when comparing 
Au+Au and Cu+Cu data from PHOBOS [i^ and is 
now extended to include Cu+Au collisions. Figure [TUI 
shows that, as in the Au+Au collisions, the Er/Nch 
ratio in the lighter colliding system is consistent with 
being independent of A^part- 

Figure [TT] shows the A^part dependence of Ebj multi¬ 
plied by r for Cu+Cu and Cu+Au collisions. Both the 
Cu+Cu data at = 200 GeV and the Cu+Cu data 

V^ivjv = 62.4 GeV increase with increasing Vpart- For 
all values of Vpart, £bj for Cu+Cu collisions at = 

200 GeV and Cu+Au at = 200 GeV are consistent 

with each other within the uncertainties of the measure¬ 
ment. With the different collision geometries taken into 
account, there is a more consistent agreement between 
the most central Cu+Cu and Cu+Au data points at 
= 200 GeV than with {dET/dr])/{0.5Npart) alone. 
Also shown for comparison are the £bj values for Au+Au 
collisions at y/s^ = 200 and 62.4 GeV, illustrating that 
Ebj is independent of the size of the system. 

VI. RESULTS FOR U+U COLLISIONS 

During the 2012 data-taking period, RBIC delivered 
U+U collisions at y/s^ =193 GeV. U+U collisions can 
provide additional information about the dynamics of the 
system [4ll - l43| by varying the collisi on g eometry of the 
nonspherical prolate uranium nuclei |45l |. Bowever, for 
this study, there is no collision geometry selection applied 
to the data. The results presented here are integrated 
over all orientations of the colliding nuclei. 

The estimate of Vpart as a function of centrality for 
U+U collisions is made using the method described pre¬ 
viously. Bowever, the U+U collisions are now modeled 
in the Glauber Monte-Carlo calculation using a deformed 
Woods Saxon distribution for the uranium nucleus to de¬ 
scribe its prolate shape, 

p(r) = po/(l + (4) 

where po is the normal nuclear density, a is the surface 
diffuseness parameter, and R' is a 0-dependent descrip¬ 
tion of the nuclear radius, 

R’ = R[1 + 132Y°{9) + /34Vi>(0)], (5) 

where V° is a Legendre Polynomial. The Woods Saxon 
parameters used are taken from a previous study with R 
= 6.81 fm, a = 0.6 fm, [32 = 0.28, and [34 = 0.093 (d^ . 
There is an additional study that presents a different set 
of parameters (i? = 6.86 fm, a = 0.42 fm, P 2 = 0.265, 
and ^4 = 0) |47l |. The two parametrizations result in 
A^part estimates that are consistent within the uncertain¬ 
ties, so the Vpart values quoted here are from the former 
parametrization 

Figure [TU] shows {dET/dr])/{0.5Nps.rt) and 
((iVch/d??)/(0.5Vpart) at midrapidity as a function 


of Vpart for U+U collisions at = 193 GeV. Also 

shown for comparison are the data for Au+Au collisions 
fo = 200 GeV. Both the U+U and the Au+Au 

data are consistent with each other for all values of 
A^part- This behavior is also observed when comparing 
Au+Au, Cu+Au, and Cu+Cu data as discussed in the 
previous section. 


VII. RESULTS FOR DEUTERON+AU AND 
^HE+AU COLLISIONS 

Measurements of dNch/dp have been publishe d by 
PBOBOS for d+Au collisions at ysJTJT = 200 GeVflTf. 
Bere those measurements are extended to include mea¬ 
surements of dEx /drj and the addition of measurements 
from ^Be+Au collisions at = 200 GeV. 

A detailed description of the method used to define 
the centrality of 200 GeV d+Au collisions using the 
PBENIX detector can be found elsewhere [d^. The same 
method was applied to define the centrality in ^Be+Au 
collisions. Figure [13] shows (dAr/dTy)/(0.5Apart) and 
(dAch/dr?)/(0.5Apart) as a function of Apart for d+Au 
and ^Be+Au collisions. Also shown are the most periph¬ 
eral Au+Au points at y/s^ = 200 GeV for compari¬ 
son. Within the uncertainties, the results for 200 GeV 
d+Au and ^Be+Au collisions are consistent with each 
other for all values of Apart- As with the heavier systems, 
the Et /Nch ratio is consistent with being independent 
of Apart within the uncertainties of the measurement as 
shown in Fig. [TH 

For minimum-bias p+p collisions at =200 GeV, 

(dAT/d??)/(0.5Apart) is 2.27 ± 0.19 GeV and 

(dAch/d77)/(0.5Apart) is 2.38 ± 0.17, where the un¬ 
certainties represent the total statistical and systematic 
uncertainties. These measurements are consistent with 
the most peripheral results from both ^Be+Au and d+Au 
collisions. The (dAch/dry)/(0.5Apart) measurement is 
also consistent with the PBOBOS measurement [i3- 


VIII. QUARK PARTICIPANT SCALING AT 
MIDRAPIDITY 

Thus far, dEx/dp and dNci^/dp have been discussed in 
terms of the dependence on the number of nucleon partic¬ 
ipants in the collision. Bere, the behavior as a function of 
the number of quark participants, A^p, will be examined. 
PBOBOS dNch/dp data for Au+Au collisions at 
= 200 and 130 GeV have been analyzed as a function of 
This analysis shows that the data at midrapidity 
are better described by scaling with A^p than with Apart 
at the top RBIC energies. A separate analysis of the 
PBOBOS dNch/dp data for Au+Au collisions extended 
down to y+Tiv = 62.4 and 19.6 GeV in terms of Nqp [i^ 
concludes that Nqp scaling better describes the data than 
Apart scaling at those lower energies. PBENIX compared 
various models of particle production and verified that 
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FIG. 12. (Color online) (di?T/rf»?)/(0.5A'part) (a) and {dNch/drj)/{0.5Npa,rt) (b) at midrapidity as a function of 
A^part for U+U collisions. Also shown are results from Au+Au collisions at = 200 GeV for comparison. The 

lines bounding the points represent the trigger efficiency uncertainty within which the points can be tilted. The 
error bars represent the remaining total statistical and systematic uncertainty. 




FIG. 13. (Color online) {dEr/drf)/{Q.bNpa,rt) (a) and {dNch/drj)/{0.5Npa,rt) (b) at midrapidity as a function of 
A^part for d+Cu and ®He+Au collisions. Also shown are results from the most peripheral Au+Au collisions at y^sTTj^T 
= 200 GeV for comparison. The lines bounding the points represent the trigger efficiency uncertainty within which 
the points can be tilted. The error bars represent the remaining total statistical and systematic uncertainty. 


Ngp scaling best describes the midrapidity (IEt / dr] mea¬ 
surements in Au+Au collisions at ^/s^ = 200 and 62.4 
GeV 0. Here these analyses are extended to include 
dET/dr] and dNch/drj measurements down to ^/s^ = 
7.7 GeV. 

The number of quark participants is estimated u sing 
a Monte Carlo Glauber model calculation method 
that has been modified to replace nucleons with con¬ 
stituent quarks Q. The nuclei are initially assembled 
by distributing the centers of the nucleons according to 


a Woods-Saxon distribution. After a nucleus is fully as¬ 
sembled, the nucleons are replaced by three quarks dis¬ 
tributed around the center of each nucleon. The quarks 
are distributed radially by sampling an empirically de¬ 
termined function: 

/(r) = r2e“'‘-2^’'(1.21466- 1.888r + 2.03r2) 

p (6) 

(1 + 1.0/r - 0.03/r^)(l + O.lSr), 

where r is the radial position of the quark in fm (4^ . The 
azimuthal position of each quark is assigned randomly to 
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FIG. 14. The Er/Nch ratio as a function of A^part for 200 GeV d+Au (a) and 200 GeV ®He+Au (b) collisions. The error bars 
represent the total statistical and systematic uncertainties. 


TABLE VI. Snmmary of the results of the linear fits to the functions dEx/dTj — asNqp +6 _b and dNch/dr) = ajvVgp +bN- 


y/^NN 

System 

as [GeV] 

bE [GeV] 

ajv 

6 jv 

200 GeV 

Au-I-Au 

0.629 ±0.021 

-6.1 ±5.4 

0.716 ±0.020 

-6.0 ±6.2 

200 GeV 

Gu-I-Au 

0.612 ±0.021 

3.4 ±2.7 

0.706 ± 0.029 

2.1 ±3.7 

200 GeV 

Cu-bCu 

0.632 ± 0.039 

1.9 ±3.9 

0.735 ± 0.040 

-1.1 ±3.9 

130 GeV 

Au-I-Au 

0.555 ±0.017 

-1.9 ±4.3 

0.635 ±0.016 

-1.6 ±4.2 

62.4 GeV 

Au-I-Au 

0.435 ±0.015 

-1.9 ±3.7 

0.499 ± 0.023 

2.2 ±5.2 

62.4 GeV 

Cu-bCu 

0.449 ± 0.026 

2.7 ± 2.8 

0.578 ± 0.043 

-0.9 ±4.5 

39 GeV 

Au-I-Au 

0.356 ±0.013 

0.8 ±3.6 

0.409 ± 0.020 

1.5 ±4.8 

27 GeV 

Au-I-Au 

0.298 ±0.010 

2.9 ±2.2 

0.357 ±0.017 

0.3 ±3.4 

19.6 GeV 

Au-I-Au 

0.264 ±0.011 

3.0 ±2.8 

0.320 ±0.016 

1.5 ±3.9 

14.5 GeV 

Au-I-Au 

0.232 ±0.010 

-1.2 ±2.5 

0.287 ±0.015 

-3.2 ±3.5 

7.7 GeV 

Au-I-Au 

0.163 ±0.007 

-1.8 ± 1.8 

0.226 ±0.017 

-2.9 ±2.9 


achieve a spherically symmetric distribution. Once all of 
the quark coordinates are determined, the center-of-mass 
of the three quark system is shifted to match the center 
position of the nucleon. The empirical function above is 
chosen such that after the center-of-mass is shifted, the 
radial distribution of the quark positions with respect to 
the nucleon center position reproduces the Fourier trans¬ 
form of the proton form factor as measured in electron- 
proton elastic scattering [50j |: 

X e-“", (7) 

where a = \/T2/rm = 4.27 fm~^ and = 0.81 fm is 
the rms charge radius of the proton 0. Once all quarks 
in both nuclei are positioned, the coordinates of the two 


^ This approach is necessary because if (^r) itself is simply 

sampled for the quark radial coordinates, the re-centering of the 
three quark system would result in a distortion of the radial 
distribution, which would then be calculated with respect to the 
center of mass of the generated system. 


nuclei are shifted relative to each other at random uni¬ 
formly in the impact parameter plane transverse to the 
beam axis. Interactions between a pair of quarks, one 
from each nucleus, occur if the distance d in this plane 
satishes the condition: 


d< 



( 8 ) 


where cr™®' is the inelastic quark-quark cross section. The 
value of 0 "“®^ is set to reproduce the known inelastic 
nucleon-nucleon cross section when running the model 
for nucleon-nucleon collisions at a given collision energy. 
The inelastic cross sections as a function of are 

taken from parametrizations of cross section measure¬ 
ments [5l| . A summary of cr™®^ as a function of is 

given in Table IIVI 

The values of midrapidity dEx/dr] and dNch/drj as 
a function of Ngp are shown in Fig. [15] for Au-|-Au 
collisions. For all collision energies, the dependence 
on Nqp is linear. When [dEr/d-q)/(Q.hNqp) and 
(dfVch/df?)/(0.57Vgp) is plotted as a function of Nqp as 
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FIG. 15. (Color online) dEr/dr] (a) and dNch/dri (b) at midrapidity as a function of Nqp for Au+Au collisions 
\/^NN ~ 200, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV. The error bars represent the total statistical and 
systematic uncertainties. 
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FIG. 16. (Color online) {dEr/dr])/{0.5Nqp) (a) and {dNch/dri)/(0.5Nqp) (b) at midrapidity as a function of Nqp 
for Au+Au collisions at = 200, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV. The lines bounding the points 

represent the trigger efficiency uncertainty within which the points can be tilted. The error bars represent the 
remaining total statistical and systematic uncertainty. 


shown in Fig. [161 the distributions are constant within 
the uncertainties of the measurement, which is not the 
case when centrality is expressed in terms of Afpart, shown 
in Fig. 121 For Au+Au collisions from = 200 to 

7.7 GeV, scaling with better describes the data than 
scaling with A^part- 

Because there is a linear dependence of dEx/dri and 
dNcY^/drj with A"qp, the data for each collision energy in 
Fig.HSlcan be fit to a straight line dEr/dr] = asNgp +bE 
and dNcYi/dr] = arrNqp +& 7 v- The extracted slopes, ar 
and Oat, represent the dEr/dr] and dNcii/dr] per quark 


participant, respectively. For all collision energies, the 
intercept of the fit at y/s//// = 0, which is kept as a 
free parameter in the fit, is consistent with zero within 
at most 1.3 standard deviations for all data sets. Fig¬ 
ure dZl shows the excitation function of the slopes for 
Au+Au collisions. The dEr/dr] data can be described 
by a second-order polynomial: ue = 0.0408 + 0.0273 x 
+ 0.0160 X The dNch/dg data 

can be described by a second-order polynomial: un = 
0.153-0.0096X Zo 5 (^/s;( 7 )+ 0.0221 x (/o5(^s);;7))2. The 
results of the linear fits for each collision energy are tab- 
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FIG. 17. (Color online) The slopes of the fit to dET/dr), ue (a), and dNch/drj, ajv (b), as a function of Nqp plotted as a function 
of for Au+Au collisions. The error bars are the uncertainties from the ht. The red line is a second-order polynomial fit 

to the data. 
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FIG. 18. (Color online) dEx/dri (a) and dNch/drj (b) at midrapidity as a function of Nqp for Cu-|-Cu and Cu+Au collisions. 
The error bars represent the total statistical and systematic uncertainties. 


ulated in Table IVllFI . 

The preference of the scaling with Nqp is also apparent 
in Cu+Cu and Cu+Au collisions. This is demonstrated 
in Fig. [T8l which shows that dEx/dri and dA^ch/dry in¬ 
creases linearly with increasing Nqp. As previously shown 
inFig.IHl (dAT/d?7)/(0.5A^part) and {dNch/drf)/[Q.bNp^^t) 


^ Note that the method of generating constituent miarks in the 
present work is slightly different than that of Ref. Q , which did 
not preserve the center of mass of the three quarks. There is a 
small effect of the different methods indicated by the small differ¬ 
ence of (di? 7 '/d» 7 } /Nqp = 0.617±0.23 GeV in Ref. Q compared 
to the present {dEj’/drj) /Nqp = 0.629 ± 0.021 GeV. 


both exhibit a distinct increase as Afpart increases for 
all three systems. This is not the case when compar¬ 
ing to Fig.[Tni which shows that {dEx/drj)/{O.bNqp) and 
{dNch/dr])/{0.5Nqp) exhibits no significant dependence 
on Nqp for all three systems. Scaling with Nqp for d+Au 
and ^He+Au collisions at = 200 GeV is shown in 

Fig. [20] along with a comparison to the most peripheral 
Au+Au collisions at = 200 GeV. As seen when 

scaled with Vpart in Fig. |T21 {dEx/dr])/(Q.^Nqp) and 
{dNch/dr])/{0.5Nqp) are also consistent with Nqp scal¬ 
ing, with the exception of {dNc^/dr])/{0.5Nqp) for d+Au 
collisions. There is no significant evidence that either 
A^part or Nqp Scaling are preferred in d+Au and ^He+Au 
collisions. 
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FIG. 19. (Color online) [dEx/drj)/ {0.5Nqp) (a) and (dNch/dri)/ (0.5Nqp) (b) at midrapidity as a function of Nqp for 
Cu+Cu and Cu+An collisions. The lines bounding the points represent the trigger efficiency uncertainty within 
which the points can be tilted. The error bars represent the remaining total statistical and systematic uncertainty. 




FIG. 20. (Color online) {dEr/dri)/ {0.5Nqp) (a) and {dNch /dr])/ {0.5Nqp) (b) at midrapidity as a function of Nqp for 
d+An and ®He+Au collisions. Shown comparison are data from the most peripheral collisions Au+An collisions 
S't \/^NN ~ 200 GeV. The lines bonnding the points represent the trigger efficiency uncertainty within which the 
points can be tilted. The error bars represent the remaining total statistical and systematic uncertainty. 


IX. SUMMARY 


Midrapidity distributions of transverse energy, 
dEx/dr], and charged particle multiplicity, dNc^/dr], 
have been measured for a variety of collision systems 
and energies, including Au+Au collisions from y/s^ = 
7.7 to 200 GeV. The centrality dependent distributions 
are presented in terms of the number of nucleon par¬ 
ticipants, Afpart! and the number of constituent quark 
participants, Ngp. The data are better described by 


scaling with Nqp than scaling with Apart- This holds for 
Au+Au collisions from = 200 GeV down to 

= 7.7 GeV, for Gu+Au collisions at = 200 GeV, 

and for Cu+Gu collisions at y/s^ = 62.4 and 200 GeV. 
Although comparisons of the data to models such as 
HIJING, parton saturation models like EKRT and KLN, 
and multiphase transport models such as AMPT are 
met with some success, a simple description using Nqp 
scaling describes the data very well. 

Some of the outstanding features of the data in¬ 
clude the following. It is observed that measurements 
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of (d£;T/c??7)/(0.5A^part) and {dNch/dri)/{0.5Nps,rt) from 
a variety of systems including Au+Au, Cu+Au, and 
Cu+Cu at = 200 GeV are all consistent with each 

other as a function of iVpart- The production of Et and 
Nch in collisions of symmetric nuclei depends only on the 
collision energy and is independent of the size of the col¬ 
liding system. The centrality dependent distributions of 
the Bjorken energy density eej show an increasing trend 
with both A^part and At = 200 GeV, eej for 

Cu-|-Au and Gu-fCu collisions are consistent with each 
other for all Vpartj again demonstrating that Et pro¬ 
duction is independent of the system size. The ratio of 
dEr/dri to dNch/dij is found to be constant as a function 
of centrality for all collision systems and energies. There 
is also only a weak dependence of this ratio as function 
of from = 7.7 to 200 GeV. Taking the ra¬ 

tio of {dET/dri)/{0.5Npart) and (dVch/(i?7)/(0.5iVpart) for 
Y^s^ = 200 GeV to 7.7 GeV shows that the shape of 
the distributions as a function of Vpart do not change 
significantly over this collision energy range. For cen¬ 
tral Au-I-Au collisions from y^F^ = 200 to 7.7 GeV, the 
value of {dET/dT])/{0.5Np^^t) and (dVch/(i77)/(0.5iVpart) 
exhibits a power law behavior as a function of y^s^^- 
Extending this observation, the Bjorken energy density 
also exhibits a power law behavior in central Au-|-Au col¬ 
lisions from Y^Sjvjv = 200 to 7.7 GeV. Also calculations of 
dEr/dr] and dNch/dr] per quark participant are observed 
to increase as yTs^ increases. 
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APPENDIX 

This Appendix contains data tables for the dET/dr] 
and dNch/drj measurements for each of the collision sys¬ 
tems. 
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TABLE VII. Transverse energy results for 200 GeV Au+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 

Centrality (A^part )_ {N^p ) dEr/drj [GeV] {dEr/dri)/{0.5Np^,t) [GeV] {dEr/dri)/ [GeV] 


0%-5% 

350.9 ± 4.7 

924.1 ± 16.2 

599.0 

± 

34.7 

3.41 

± 

0.20 

1.30 

± 0.08 

5%-10% 

297.0 ± 6.6 

782.6 ± 15.3 

498.7 

± 

28.9 

3.30 

± 

0.21 

1.25 

± 0.08 

10%-15% 

251.0 ± 7.3 

644.6 ± 14.5 

403.0 

± 

25.0 

3.21 

± 

0.22 

1.25 

± 0.08 

15%-20% 

211.0 ± 7.3 

532.9 ± 12.3 

332.5 

± 

21.2 

3.15 

± 

0.23 

1.25 

± 0.08 

20%-25% 

176.3 ± 7.0 

437.5 ± 10.4 

273.6 

± 

18.6 

3.10 

± 

0.24 

1.25 

± 0.09 

25%-30% 

146.8 ± 7.1 

356.8 ± 12.2 

223.4 

± 

16.4 

3.04 

± 

0.27 

1.25 

± 0.10 

30%-35% 

120.9 ± 7.0 

288.3 ± 11.0 

180.8 

± 

14.3 

2.99 

± 

0.29 

1.25 

± 0.11 

35%-40% 

98.3 ± 6.8 

229.7 ± 9.2 

144.5 

± 

12.6 

2.94 

± 

0.33 

1.26 

± 0.12 

40%-45% 

78.7 ± 6.1 

181.0 ± 6.8 

113.9 

± 

10.9 

2.89 

± 

0.36 

1.26 

± 0.13 

45%-50% 

61.9 ± 5.2 

141.1 ± 5.3 

88.3 

± 

9.3 

2.85 

± 

0.38 

1.25 

± 0.14 

50%-55% 

47.6 ± 4.9 

107.6 ± 5.5 

67.1 

± 

8.1 

2.82 

± 

0.45 

1.25 

± 0.16 

55%-60% 

35.6 ± 5.1 

77.5 ± 6.8 

50.0 

± 

6.7 

2.81 

± 

0.55 

1.29 

± 0.21 


TABLE VIII. Charged particle multiplicity results for 200 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(A^part ) 

(A,p ) 

dNch/dri 

(dAch/dr?)/(0.5Apart) (dNch/drj)/(0.5Nqp) 

0%-5% 

350.9 ± 4.7 

924.1 ± 16.2 

687.4 ± 36.6 


3.92 ± 0.22 

1.49 

O 

o 

00 

5%-10% 

297.9 ± 6.6 

782.6 ± 15.3 

560.4 ± 27.9 


3.77 ± 0.21 

1.43 

± 0.08 

10%-15% 

251.0 ± 7.3 

644.6 ± 14.5 

456.8 ± 22.3 


3.64 ± 0.21 

1.42 

± 0.08 

15%-20% 

211.0 ± 7.3 

532.9 ± 12.3 

371.5 ± 18.2 


3.52 ± 0.21 

1.39 

± 0.08 

20%-25% 

176.3 ± 7.0 

437.5 ± 10.4 

302.5 ± 15.8 


3.43 ± 0.22 

1.38 

± 0.08 

25%-30% 

146.8 ± 7.1 

356.8 ± 12.2 

245.6 ± 13.8 


3.35 ± 0.25 

1.38 

± 0.09 

30%-35% 

120.9 ± 7.0 

288.3 ± 11.0 

197.2 ± 12.2 


3.26 ± 0.28 

1.37 

± 0.10 

35%-40% 

98.3 ± 6.8 

229.7 ± 9.2 

156.4 ± 10.9 


3.18 ± 0.31 

1.36 

± 0.11 

40%-45% 

78.7 ± 6.1 

181.0 ± 6.8 

123.5 ± 9.6 


3.14 ± 0.34 

1.36 

± 0.12 

45%-50% 

61.9 ± 5.2 

141.1 ± 5.3 

95.3 ± 8.6 


3.08 ± 0.38 

1.35 

± 0.13 

50%-55% 

47.6 ± 4.9 

107.6 ± 5.5 

70.9 ± 7.6 


2.98 ± 0.44 

1.32 

± 0.16 

55%-60% 

35.6 ± 5.1 

77.5 ± 6.8 

52.2 ± 6.5 


2.93 ± 0.56 

1.35 

± 0.20 

TABLE IX. 

Transverse energy 

results for 130 GeV Au+Au collisions. 

The uncertainties include the total statistical and 

systematic uncertainties. 







Centrality 

(Apart ) 

(A,p ) dET/d-n [GeV] (dAT/dr?)/(0.5Apart) [GeV] 

{dET/d'n)/{0.5Nqp) [GeV] 

0%-5% 

347.7 ± 10.0 

914.1 ± 22.6 

522.8 ± 27.3 

3.01 

± 0.18 

1.14 ± 

0.07 

5%-10% 

294.0 ± 8.9 

773.3 ± 20.3 

425.2 ± 22.5 

2.89 

± 0.18 

1.10 ± 

0.07 

10%-15% 

249.5 ± 8.0 

633.4 ± 19.4 

349.0 ± 19.0 

2.80 

± 0.18 

1.10 ± 

0.07 

15%-20% 

211.0 ± 7.2 

522.6 ± 18.3 

287.2 ± 16.5 

2.72 

± 0.18 

1.10 ± 

0.07 

20%-25% 

178.6 ± 6.6 

431.5 ± 19.0 

237.1 ± 14.5 

2.66 

± 0.19 

1.10 ± 

0.08 

25%-30% 

149.7 ± 6.0 

353.3 ± 15.9 

191.3 ± 12.5 

2.56 

± 0.20 

1.08 ± 

0.09 

30%-35% 

124.8 ± 5.5 

283.0 ± 13.2 

153.9 ± 11.0 

2.47 

± 0.21 

1.09 ± 

0.09 

35%-40% 

102.9 ± 5.1 

225.3 ± 11.0 

121.8 ± 9.4 

2.37 

± 0.22 

1.08 ± 

0.10 

40%-45% 

83.2 ± 4.7 

179.1 ± 8.8 

96.0 ± 8.8 

2.31 

± 0.25 

1.07 ± 

0.11 

45%-50% 

66.3 ± 4.3 

137.1 ± 7.1 

73.3 ± 7.3 

2.21 

± 0.26 

1.07 ± 

0.12 

50%-55% 

52.1 ± 4.0 

101.6 ± 6.5 

55.5 ± 6.5 

2.13 

± 0.30 

1.09 ± 

0.15 

55%-60% 

40.1 ± 3.8 

74.6 ± 7.3 

41.0 ± 5.5 

2.04 

± 0.34 

1.10 ± 

0.18 
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TABLE X. Charged particle multiplicity results for 130 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(Apart ) 


(A,p ) 

dNch/dr] 

(dAch/d??)/(0.5Apart) 

{dNch/dr])/{0.5Nqp) 

0%-5% 

347.7 ± 10.0 


914.1 ± 22.6 601.8 ± 28.4 


3.46 ± 

0.19 

1.32 

± 0.07 

5%-10% 

294.0 ± 8.9 


773.3 ± 20.3 488.5 ± 21.6 


3.32 ± 

0.18 

1.26 

± 0.07 

10%-15% 

249.5 ± 8.0 


633.4 ± 19.4 402.7 ± 17.4 


3.23 ± 

0.17 

1.27 

± 0.07 

15%-20% 

211.0 ± 7.2 


522.6 ± 18.3 328.8 ± 15.2 


3.12 ± 

0.18 

1.26 

± 0.07 

20%-25% 

178.6 ± 6.6 


431.5 ± 19.0 270.5 ± 12.8 


3.03 ± 

0.18 

1.25 

± 0.08 

25%-30% 

149.7 ± 6.0 


353.3 ± 15.9 219.3 ± 11.4 


2.93 ± 

0.19 

1.24 

± 0.09 

30%-35% 

124.8 ± 5.5 


283.0 ± 13.2 175.7 ± 10.3 


2.82 ± 

0.21 

1.24 

± 0.09 

35%-40% 

102.9 ± 5.1 


225.3 ± 11.0 139.0 ± 9.1 


2.70 ± 

0.22 

1.23 

± 0.10 

40%-45% 

83.2 ± 4.7 


179.1 ± 8.8 

1 109.4 ± 8.4 


2.63 ± 

0.25 

1.22 

± 0.11 

45%-50% 

66.3 ± 4.3 


137.1 ± 7.1 

84.1 ± 7.0 


2.54 ± 

0.27 

1.23 

± 0.12 

50%-55% 

52.1 ± 4 


101.6 ± 6.5 

I 64.3 ± 6.3 


2.47 ± 

0.31 

1.27 

± 0.15 

55%-60% 

40.1 ± 3.8 


74.6 ± 7.3 

48.4 ± 5.4 


2.41 ± 

0.35 

1.30 

± 0.19 

TABLE XL 

Transverse energy results for 62.4 GeV Au+Au collisions. 

The uncertainties include the total statistical and 

systematic uncertainties. 









Centrality 

(A^part ) 

(A^p ) dEr/dTj [GeV] (dAT/d??)/(0.5Apart) [GeV] 

(dAT/dr?)/(0.5A,p) [GeV] 

0%-5% 

342.6 ± 4.9 

891.7 

± 26.6 

389.7 ± 23.5 

2.27 

± 0.14 


0.87 ± 

0.06 

5%-10% 

291.3 ± 7.3 

730.7 

± 24.1 

320.5 ± 19.3 

2.20 

± 0.14 


0.88 ± 

0.06 

10%-15% 

244.5 ± 8.9 

600.6 

± 21.5 

260.6 ± 15.7 

2.13 

± 0.15 


0.87 ± 

0.06 

15%-20% 

205.0 ± 9.6 

493.4 

± 19.6 

212.1 ± 12.8 

2.07 

± 0.16 


0.86 ± 

0.06 

20%-25% 

171.3 ± 8.9 

403.8 

± 18.5 

171.9 ± 10.4 

2.01 

± 0.16 


0.85 ± 

0.06 

25%-30% 

142.2 ± 8.5 

327.0 

± 16.7 

138.6 ± 8.36 

1.95 

± 0.17 


0.85 ± 

0.07 

30%-35% 

116.7 ± 8.9 

261.7 

± 15.7 

110.4 ± 6.67 

1.90 

± 0.18 


0.85 ± 

0.07 

35%-40% 

95.2 ± 7.7 

206.9 

± 14.3 

86.9 ± 5.25 

1.83 

± 0.19 


0.84 ± 

0.08 

40%-45% 

76.1 ± 7.7 

161.4 

± 13.3 

67.3 ± 4.08 

1.78 

± 0.21 


0.84 ± 

0.09 

45%-50% 

59.9 ± 6.9 

123.5 

± 13.2 

51.2 ± 3.12 

1.73 

± 0.22 


0.84 ± 

0.10 

50%-55% 

46.8 ± 5.2 

92.4 ± 11.2 

38.4 ± 2.33 

1.65 

± 0.21 


0.84 ± 

0.11 

55%-60% 

35.8 ± 4.6 

67.8 ± 9.0 

28.5 ± 1.72 

1.59 

± 0.23 


0.84 ± 

0.12 


TABLE XII. Charged particle multiplicity results for 62.4 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(Apart ) 

(A,p 

) 

dAch/d?? 

(dAch/dr;) / (0.5 Apart) 

(dAch/dr;)/(0.5A,p) 

0%-5% 

342.6 ± 4.9 

891.7 ± 

26.6 

447.5 ± 38.9 

2.61 

± 

0.23 

1.00 

± 

0.09 

5%-10% 

291.3 ± 7.3 

730.7 ± 

24.1 

367.4 ± 31.6 

2.52 

± 

0.23 

1.01 

± 

0.09 

10%-15% 

244.5 ± 8.9 

600.6 ± 

21.5 

301.8 ± 25.8 

2.47 

± 

0.23 

1.01 

± 

0.09 

15%-20% 

205.0 ± 9.6 

493.4 ± 

19.6 

248.0 ± 21.0 

2.42 

± 

0.23 

1.01 

± 

0.09 

20%-25% 

171.3 ± 8.9 

403.8 ± 

18.5 

203.0 ± 17.1 

2.37 

± 

0.24 

1.01 

± 

0.10 

25%-30% 

142.2 ± 8.5 

327.0 ± 

16.7 

165.1 ± 13.8 

2.32 

± 

0.24 

1.01 

± 

0.10 

30%-35% 

116.7 ± 8.9 

261.7 ± 

15.7 

133.0 ± 11.1 

2.28 

± 

0.26 

1.02 

± 

0.10 

35%-40% 

95.2 ± 7.7 

206.9 ± 

14.3 

105.9 ± 8.76 

2.22 

± 

0.26 

1.02 

± 

0.11 

40%-45% 

76.1 ± 7.7 

161.4 ± 

13.3 

83.0 ± 6.83 

2.18 

± 

0.28 

1.03 

± 

0.12 

45%-50% 

59.9 ± 6.9 

123.5 ± 

13.2 

63.9 ± 5.24 

2.13 

± 

0.30 

1.03 

± 

0.14 

50%-55% 

46.8 ± 5.2 

92.4 ± 11.2 

48.4 ± 3.95 

2.07 

± 

0.29 

1.05 

± 

0.15 

55%-60% 

35.8 ± 4.6 

67.8 ± 9.0 

35.8 ± 2.92 

2.00 

± 

0.30 

1.06 

± 

0.16 
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TABLE XIII. Transverse energy results for 39 GeV Au+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 


Centrality 

(Apart ) 

{Nqp 

) 

dEr/dri [GeV] 

(dAT/dr?)/(0.5Apart) [GeV] 

{dET/d'n)/{0.5Ngp) [GeV] 

0%-5% 

340.0 ± 7.4 

874.6 

± 

42.0 

303.8 ± 18.2 

1.79 

± 

0.11 

0.69 

± 

0.05 

5%-10% 

289.6 ± 8.1 

726.7 

± 

36.7 

262.1 ± 15.7 

1.81 

± 

0.12 

0.72 

± 

0.06 

10%-15% 

244.1 ± 6.4 

599.1 

± 

26.8 

216.6 ± 13.0 

1.77 

± 

0.12 

0.72 

± 

0.05 

15%-20% 

206.5 ± 6.3 

496.9 

± 

23.7 

178.5 ± 10.7 

1.73 

± 

0.12 

0.72 

± 

0.06 

20%-25% 

174.1 ± 6.3 

410.4 

± 

20.9 

146.9 ± 8.8 

1.69 

± 

0.12 

0.72 

± 

0.06 

25%-30% 

145.8 ± 6.2 

336.8 

± 

22.2 

120.4 ± 7.2 

1.65 

± 

0.12 

0.72 

± 

0.06 

30%-35% 

120.8 ± 7.5 

273.0 

± 

18.1 

97.7 ± 5.8 

1.62 

± 

0.14 

0.72 

± 

0.06 

35%-40% 

98.6 ± 6.4 

217.6 

± 

15.1 

78.5 ± 4.7 

1.59 

± 

0.14 

0.72 

± 

0.07 

40%-45% 

79.8 ± 6.0 

172.0 

± 

14.1 

62.3 ± 3.7 

1.56 

± 

0.15 

0.72 

± 

0.07 

45%-50% 

63.9 ± 5.8 

134.3 

± 

13.1 

48.6 ± 2.9 

1.52 

± 

0.17 

0.72 

± 

0.08 

50%-55% 

50.3 ± 5.5 

103.1 

± 

13.5 

37.3 ± 2.2 

1.48 

± 

0.18 

0.72 

± 

0.10 


TABLE XIV. Charged particle multiplicity results for 39 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(Apart ) 

(A,p 

) 

dAch/d?? 

(dAch/d?))/(0.5 Apart ) 

(dAch/dr;)/(0.5A,p) 

0%-5% 

340.0 ± 7.4 

874.6 

± 

42.0 

363.2 ± 31.6 

2.14 

± 

0.19 

0.83 

± 

0.08 

5%-10% 

289.6 ± 8.1 

726.7 

± 

36.7 

297.8 ± 25.8 

2.06 

± 

0.19 

0.82 

± 

0.08 

10%-15% 

244.1 ± 6.4 

599.1 

± 

26.8 

246.6 ± 21.3 

2.02 

± 

0.18 

0.82 

± 

0.08 

15%-20% 

206.5 ± 6.3 

496.9 

± 

23.7 

204.4 ± 17.5 

1.98 

± 

0.18 

0.82 

± 

0.08 

20%-25% 

174.1 ± 6.3 

410.4 

± 

20.9 

168.9 ± 14.4 

1.94 

± 

0.18 

0.82 

± 

0.08 

25%-30% 

145.8 ± 6.2 

336.8 

± 

22.2 

138.3 ± 11.8 

1.90 

± 

0.18 

0.82 

± 

0.09 

30%-35% 

120.8 ± 7.5 

273.0 

± 

18.1 

112.6 ± 9.6 

1.86 

± 

0.20 

0.83 

± 

0.09 

35%-40% 

98.6 ± 6.4 

217.6 

± 

15.1 

90.6 ± 7.7 

1.84 

± 

0.20 

0.83 

± 

0.09 

40%-45% 

79.8 ± 6.0 

172.0 

± 

14.1 

72.1 ± 6.1 

1.81 

± 

0.20 

0.84 

± 

0.10 

45%-50% 

63.9 ± 5.8 

134.3 

± 

13.1 

56.8 ± 4.8 

1.78 

± 

0.22 

0.85 

± 

0.11 

50%-55% 

50.3 ± 5.5 

103.1 

± 

13.5 

43.7 ± 3.7 

1.73 

± 

0.24 

0.85 

± 

0.13 


TABLE XV. Transverse energy results for 27 GeV Au+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 


Centrality 

(A^part ) 

{Ngp 

) 

dEr/dri [GeV] 

(dAT/dr?)/(0.5Apart) [GeV] 

(dAT/dr?)/(0.5A,p) [GeV] 

0%-5% 

338.9 ± 3.1 

863.7 

± 

23.5 

265.6 ± 16.4 

1.57 

± 

0.10 

0.62 

± 

0.04 

5%-10% 

288.8 ± 4.7 

718.8 

± 

22.7 

217.3 ± 13.4 

1.50 

± 

0.10 

0.61 

± 

0.04 

10%-15% 

244.3 ± 6.5 

595.0 

± 

23.7 

179.7 ± 11.1 

1.47 

± 

0.10 

0.60 

± 

0.04 

15%-20% 

205.7 ± 5.8 

490.7 

± 

19.4 

148.9 ± 9.2 

1.45 

± 

0.10 

0.61 

± 

0.04 

20%-25% 

173.0 ± 5.5 

404.6 

± 

16.7 

122.8 ± 7.6 

1.42 

± 

0.10 

0.61 

± 

0.05 

25%-30% 

144.6 ± 6.2 

330.8 

± 

17.7 

100.7 ± 6.2 

1.39 

± 

0.10 

0.61 

± 

0.05 

30%-35% 

119.4 ± 6.1 

267.4 

± 

16.2 

81.9 ± 5.1 

1.37 

± 

0.11 

0.61 

± 

0.05 

35%-40% 

97.6 ± 5.8 

213.6 

± 

14.6 

65.8 ± 4.1 

1.35 

± 

0.12 

0.62 

± 

0.06 

40%-45% 

77.9 ± 5.7 

166.0 

± 

13.7 

52.1 ± 3.2 

1.34 

± 

0.13 

0.63 

± 

0.06 

45%-50% 

60.8 ± 6.0 

125.9 

± 

13.8 

40.8 ± 2.5 

1.34 

± 

0.16 

0.65 

± 

0.08 


TABLE XVI. Charged particle multiplicity results for 27 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 

Centrality_(Aipart )_ {Ngp ) _ dN^h/dr] _(cLVch/d?7)/(0.5jVpart)_ {dN^h/drj)/{Q.^Ngp) 


0%-5% 

338.9 

± 

3.1 

863.7 

± 

23.5 

321.2 

± 

28.1 

1.90 

± 

0.17 

0.74 

± 

0.07 

5%-10% 

288.8 

± 

4.7 

718.8 

± 

22.7 

258.7 

± 

22.5 

1.79 

± 

0.16 

0.72 

± 

0.07 

10%-15% 

244.3 

± 

6.5 

595.0 

± 

23.7 

212.6 

± 

18.5 

1.74 

± 

0.16 

0.72 

± 

0.07 

15%-20% 

205.7 

± 

5.8 

490.7 

± 

19.4 

175.0 

± 

15.1 

1.70 

± 

0.16 

0.71 

± 

0.07 

20%-25% 

173.0 

± 

5.5 

404.6 

± 

16.7 

143.5 

± 

12.4 

1.66 

± 

0.15 

0.71 

± 

0.07 

25%-30% 

144.6 

± 

6.2 

330.8 

± 

17.7 

116.7 

± 

10.0 

1.61 

± 

0.16 

0.71 

± 

0.07 

30%-35% 

119.4 

± 

6.1 

267.4 

± 

16.2 

94.2 

± 

8.1 

1.58 

± 

0.16 

0.70 

± 

0.07 

35%-40% 

97.6 

± 

5.8 

213.6 

± 

14.6 

75.0 

± 

6.4 

1.54 

± 

0.16 

0.70 

± 

0.08 
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TABLE XVII. Transverse energy results for 19.6 GeV Au+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 


Centrality 

(Apart ) 

(A(jp 

) 

dEr/dT) [GeV] 

(dAT/d??)/(0.5Apart) [GeV] 

(dAT/d?7)/(0.5A,p) [GeV] 

0%-5% 

338.5 ± 4.4 

858.8 

± 

27.7 

233.1 ± 15.3 

1.38 

± 

0.09 

0.54 

± 

0.04 

5%-10% 

288.3 ± 6.0 

714.1 

± 

25.0 

190.7 ± 12.5 

1.32 

± 

0.09 

0.53 

± 

0.04 

10%-15% 

242.4 ± 6.1 

587.3 

± 

23.2 

157.8 ± 10.3 

1.30 

± 

0.09 

0.54 

± 

0.04 

15%-20% 

204.3 ± 5.7 

485.0 

± 

20.0 

130.8 ± 8.6 

1.28 

± 

0.09 

0.54 

± 

0.04 

20%-25% 

172.4 ± 7.3 

401.4 

± 

22.0 

108.2 ± 7.1 

1.25 

± 

0.10 

0.54 

± 

0.05 

25%-30% 

143.5 ± 6.6 

326.9 

± 

18.4 

88.9 ± 5.8 

1.24 

± 

0.10 

0.54 

± 

0.05 

30%-35% 

117.9 ± 6.7 

262.5 

± 

17.7 

72.5 ± 4.8 

1.23 

± 

0.11 

0.55 

± 

0.05 

35%-40% 

95.7 ± 6.9 

208.0 

± 

17.0 

58.5 ± 3.8 

1.22 

± 

0.12 

0.56 

± 

0.06 

40%-45% 

77.4 ± 5.7 

164.2 

± 

13.5 

46.6 ± 3.1 

1.21 

± 

0.12 

0.57 

± 

0.06 

45%-50% 

61.7 ± 4.8 

127.7 

± 

11.3 

36.6 ± 2.4 

1.19 

± 

0.12 

0.57 

± 

0.06 


TABLE XVIII. Charged particle multiplicity results for 19.6 GeV Au+Au collisions. The uncertainties include the total 
statistical and systematic uncertainties. 


Centrality 

(Apart ) 



(A,p 

) 

dAch/d?7 

(dAch/dr?)/(0.5Apart 

) (dAah/d77)/(0.5A,p) 

0%-5% 

338.5 ± 4.4 


858.8 ± 

27.7 

285.3 ± 25.1 


1.69 ± 0.15 

0.66 

± 0.06 

5%-10% 

288.3 ± 6.0 


714.1 ± 

25.0 

229.3 ± 20.1 


1.59 ± 0.14 

0.64 

± 0.06 

10%-15% 

242.4 ± 6.1 


587.3 ± 

23.2 

188.8 ± 16.5 


1.56 ± 0.14 

0.64 

± 0.06 

15%-20% 

204.3 ± 5.7 


485.0 ± 

20.0 

155.7 ± 13.5 


1.52 ± 0.14 

0.64 

± 0.06 

20%-25% 

172.4 ± 7.3 


401.4 ± 

22.0 

128.2 ± 11.1 


1.49 ± 0.14 

0.64 

± 0.07 

25%-30% 

143.5 ± 6.6 


326.9 ± 

18.4 

104.8 ± 9.1 


1.46 ± 0.14 

0.64 

± 0.07 

30%-35% 

117.9 ± 6.7 


262.5 ± 

17.7 

85.1 ± 7.4 


1.44 ± 0.15 

0.65 

± 0.07 

35%-40% 

95.7 ± 6.9 


208.0 ± 

17.0 

68.4 ± 5.9 


1.43 ± 0.16 

0.66 

± 0.08 

40%-45% 

77.4 ± 5.7 


164.2 ± 

13.5 

54.3 ± 4.7 


1.40 ± 0.16 

0.66 

± 0.08 

45%-50% 

61.7 ± 4.8 


127.7 ± 

11.3 

42.4 ± 3.7 


1.37 ± 0.16 

0.66 

± 0.08 

TABLE XIX. 

Transverse energy results for 

14.5 GeV Au+Au collisions. 

The uncertainties include the total statistical and 

systematic uncertainties. 










Centrality 

(Apart ) 

(A(jp 

) 

dEr/dr, [GeV] (dAT/dr?)/(0.5Apart) [GeV] 

(dAT/d?7)/(0.5A,p) [GeV] 

0%-5% 

337.3 ± 4.2 

852.0 

± 

27.5 

200.4 ± 14.0 

1.19 

± 0.08 

0.47 ± 

0.04 

5%-10% 

287.7 ± 4.9 

710.1 

± 

23.4 


164.0 ± 11.5 

1.14 

± 0.08 

0.46 ± 

0.04 

10%-15% 

242.5 ± 5.5 

585.6 

± 

22.0 


134.9 ± 9.4 

1.11 

± 0.08 

0.46 ± 

0.04 

15%-20% 

205.1 ± 5.9 

485.5 

± 

19.7 


111.0 ± 7.8 

1.08 

± 0.08 

0.46 ± 

0.04 

20%-25% 

172.6 ± 6.4 

400.4 

± 

19.6 


91.1 ± 6.4 

1.06 

± 0.08 

0.46 ± 

0.04 

25%-30% 

143.6 ± 7.8 

325.9 

± 

21.7 


74.4 ± 5.2 

1.04 

± 0.09 

0.46 ± 

0.04 

30%-35% 

119.2 ± 7.2 

264.9 

± 

19.2 


60.2 ± 4.2 

1.01 

± 0.09 

0.45 ± 

0.05 

35%-40% 

98.3 ± 5.8 

213.7 

± 

14.8 


48.2 ± 3.4 

0.98 

± 0.09 

0.45 ± 

0.04 

40%-45% 

80.2 ± 5.6 

170.2 

± 

13.6 


38.2 ± 2.7 

0.95 

± 0.09 

0.45 ± 

0.05 

45%-50% 

63.9 ± 4.7 

132.2 

± 

11.0 


29.7 ± 2.1 

0.93 

± 0.09 

0.45 ± 

0.05 


TABLE XX. Gharged particle multiplicity results for 14.5 GeV Au+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(A^part ) 

(A,p 

) 

dAch/d?7 

(dAch/dr?) / (0.5 Apart) 

(dAch/d77)/(0.5A5p) 

0%-5% 

337.3 ± 4.2 

852.0 

± 

27.5 

250.9 ± 22.2 

1.49 

± 

0.13 

0.59 

± 

0.06 

5%-10% 

287.7 ± 4.9 

710.1 

± 

23.4 

201.2 ± 17.7 

1.40 

± 

0.13 

0.57 

± 

0.05 

10%-15% 

242.5 ± 5.5 

585.6 

± 

22.0 

164.5 ± 14.5 

1.36 

± 

0.12 

0.56 

± 

0.05 

15%-20% 

205.1 ± 5.9 

485.5 

± 

19.7 

134.7 ± 11.8 

1.31 

± 

0.12 

0.56 

± 

0.05 

20%-25% 

172.6 ± 6.4 

400.4 

± 

19.6 

110.0 ± 9.6 

1.28 

± 

0.12 

0.55 

± 

0.06 

25%-30% 

143.6 ± 7.8 

325.9 

± 

21.7 

89.4 ± 7.8 

1.25 

± 

0.13 

0.55 

± 

0.06 

30%-35% 

119.2 ± 7.2 

264.9 

± 

19.2 

72.0 ± 6.3 

1.21 

± 

0.13 

0.54 

± 

0.06 

35%-40% 

98.3 ± 5.8 

213.7 

± 

14.8 

57.4 ± 5.0 

1.17 

± 

0.12 

0.55 

± 

0.06 

40%-45% 

80.2 ± 5.6 

170.2 

± 

13.6 

45.2 ± 3.9 

1.13 

± 

0.13 

0.53 

± 

0.06 

45%-50% 

63.9 ± 4.7 

132.2 

± 

11.0 

34.9 ± 3.0 

1.09 

± 

0.12 

0.53 

± 

0.06 




















27 


TABLE XXL Transverse energy results for 
systematic uncertainties. 

7.7 GeV Au+Au collisions. 

The uncertainties include the total statistical and 

Centrality 

(Apart ) 

(A,p) 

dEr/dri [GeV] {dET/d'n)/{0.5Np^rt) 

[GeV] 

(dAT/d?7)/(0.5A,p) [GeV] 

0%-5% 

332.1 ± 5.4 

830.4 ± 33.9 

137.7 ± 9.1 

0.83 ± 0.06 


0.33 ± 0.03 

5%-10% 

283.2 ± 5.9 

692.3 ± 27.0 

114.3 ± 7.5 

0.81 ± 0.06 


0.33 ± 0.03 

10%-15% 

240.1 ± 5.7 

574.4 ± 24.0 

93.3 ± 6.2 

0.78 ± 0.05 


0.33 ± 0.03 

15%-20% 

204.1 ± 5.7 

479.0 ± 20.6 

76.2 ± 5.0 

0.75 ± 0.05 


0.32 ± 0.03 

20%-25% 

172.9 ± 6.7 

398.0 ± 19.8 

62.0 ± 4.1 

0.72 ± 0.05 


0.31 ± 0.03 

25%-30% 

145.5 ± 7.2 

328.1 ± 19.8 

50.0 ± 3.3 

0.69 ± 0.06 


0.30 ± 0.03 

30%-35% 

121.0 ± 7.3 

267.1 ± 19.0 

40.1 ± 2.6 

0.66 ± 0.06 


0.30 ± 0.03 

35%-40% 

98.2 ± 7.0 

211.6 ± 17.8 

31.8 ± 2.1 

0.66 ± 0.06 


0.30 ± 0.03 

40%-45% 

78.8 ± 6.7 

165.6 ± 16.3 

24.8 ± 1.6 

0.63 ± 0.07 


0.30 ± 0.04 

45%-50% 

61.8 ± 6.5 

126.4 ± 14.7 

19.2 ± 1.3 

0.62 ± 0.08 


0.30 ± 0.04 

TABLE XXII. Charged particle multiplicity results for 7.7 GeV Au-|-Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 

Centrality 

(Apart ) 

(A,p 

) dAch/d?? 

(dAch/d?;) / (0.5 Apart ) 

(dAch/d77)/(0.5A5p) 

0%-5% 

332.1 ± 5.4 

830.4 ± 

33.9 192.4 ± 16.9 

1.16 ± 

0.10 

0.46 ± 0.04 

5%-10% 

283.2 ± 5.9 

692.3 ± 

27.0 159.2 ± 14.0 

1.12 ± 

0.10 

0.46 ± 0.04 

10%-15% 

240.1 ± 5.7 

574.4 ± 

24.0 129.3 ± 11.3 

1.08 ± 

0.10 

0.45 ± 0.04 

15%-20% 

204.1 ± 5.7 

479.0 ± 

20.6 105.4 ± 9.2 

1.03 ± 

0.09 

0.44 ± 0.04 

20%-25% 

172.9 ± 6.7 

398.0 ± 

19.8 85.6 ± 7.5 

0.99 ± 

0.09 

0.43 ± 0.04 

25%-30% 

145.5 ± 7.2 

328.1 ± 

19.8 68.8 ± 6.0 

0.95 ± 

0.09 

0.42 ± 0.04 

30%-35% 

121.0 ± 7.3 

267.1 ± 

19.0 55.0 ± 4.8 

0.91 ± 

0.10 

0.41 ± 0.05 

35%-40% 

98.2 ± 7.0 

211.6 ± 

17.8 43.5 ± 3.8 

0.89 ± 

0.10 

0.41 ± 0.05 

40%-45% 

78.8 ± 6.7 

165.6 ± 

16.3 33.9 ± 3.0 

0.86 ± 

0.11 

0.41 ± 0.05 

45%-50% 

61.8 ± 6.5 

126.4 ± 

14.7 26.1 ± 2.3 

0.85 ± 

0.11 

0.41 ± 0.06 

TABLE XXIII. Transverse energy results for 200 GeV Gu-|-Gu collisions, 
systematic uncertainties. 

The uncertainties include the total statistical and 

Gentrality 

(Apart ) 

(A,p) 

dEr/dri [GeV] {dET/dr,)/{0.5Np^rt) 

[GeV] 

(dAT/d?7)/(0.5A,p) [GeV] 

0%-5% 

105.6 ± 2.5 

254.3 ± 11.8 

166.8 ± 13.2 

3.16 ± 0.26 


1.31 ± 0.12 

5%-10% 

93.1 ± 3.0 

219.0 ± 11.4 

139.9 ± 11.1 

3.01 ± 0.26 


1.28 ± 0.12 

10%-15% 

80.1 ± 2.4 

183.6 ± 8.6 

117.1 ± 9.3 

2.92 ± 0.25 


1.28 ± 0.12 

15%-20% 

68.4 ± 2.5 

153.0 ± 7.7 

97.9 ± 7.8 

2.86 ± 0.25 


1.28 ± 0.12 

20%-25% 

58.4 ± 2.3 

127.7 ± 7.0 

81.6 ± 6.5 

2.80 ± 0.25 


1.28 ± 0.12 

25%-30% 

49.2 ± 2.1 

104.9 ± 5.7 

67.8 ± 5.4 

2.76 ± 0.25 


1.29 ± 0.12 

30%-35% 

41.3 ± 2.2 

86.0 ± 5.8 

56.1 ± 4.4 

2.72 ± 0.26 


1.30 ± 0.13 

35%-40% 

34.3 ± 2.0 

69.8 ± 5.0 

46.0 ± 3.6 

2.68 ± 0.26 


1.32 ± 0.14 

40%-45% 

28.1 ± 1.8 

55.9 ± 4.3 

37.5 ± 3.0 

2.67 ± 0.27 


1.34 ± 0.15 


TABLE XXIV. Charged particle multiplicity results for 200 GeV Cu+Cu collisions. The uncertainties include the total statis¬ 
tical and systematic uncertainties. 


Gentrality 

(Apart ) 

(A,p ) 

dNch/dri 

(dAch/d?))/(0.5 Apart ) 

(dAch/dr;)/ {0.5Nqp) 

0%-5% 

105.6 ± 2.5 

254.3 ± 11.8 

192.6 ± 13.9 

3.65 ± 0.28 

1.51 ± 0.13 

5%-10% 

93.1 ± 3.0 

219.0 ± 11.4 

160.1 ± 11.5 

3.44 ± 0.27 

1.46 ± 0.13 

10%-15% 

80.1 ± 2.4 

183.6 ± 8.6 

132.8 ± 9.5 

3.32 ± 0.26 

1.45 ± 0.12 

15%-20% 

68.4 ± 2.5 

153.0 ± 7.7 

110.2 ± 7.9 

3.22 ± 0.26 

1.44 ± 0.12 

20%-25% 

58.4 ± 2.3 

127.7 ± 7.0 

91.3 ± 6.5 

3.13 ± 0.25 

1.43 ± 0.13 

25%-30% 

49.2 ± 2.1 

104.9 ± 5.7 

75.2 ± 5.3 

3.06 ± 0.25 

1.43 ± 0.13 

30%-35% 

41.3 ± 2.2 

86.0 ± 5.8 

61.7 ± 4.4 

2.99 ± 0.27 

1.43 ± 0.14 

35%-40% 

34.3 ± 2.0 

69.8 ± 5.0 

50.2 ± 3.5 

2.93 ± 0.27 

1.44 ± 0.14 

40%-45% 

28.1 ± 1.8 

55.9 ± 4.3 

40.6 ± 2.9 

2.89 ± 0.28 

1.45 ± 0.15 
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TABLE XXV. Transverse energy results for 62.4 GeV Cu+Cu collisions. The uncertainties include the total statistical and 
systematic uncertainties. 

Centrality (Vpart )_ {N^p ) dEr/drj [GeV] {dEr/drj)/{0.5Np^,t) [GeV] {dET/dr^)/{0.5Ngp) [GeV] 


0%-5% 

100.5 

± 

4.5 

229.3 

± 

8.5 

107.6 

± 

6.5 

2.14 

± 

0.16 

0.94 

± 

0.07 

5%-10% 

88.3 

± 

4.8 

197.8 

± 

15.0 

93.6 

± 

5.6 

2.12 

± 

0.17 

0.95 

± 

0.09 

10%-15% 

78.2 

± 

4.3 

171.7 

± 

25.2 

79.3 

± 

4.8 

2.03 

± 

0.17 

0.92 

± 

0.15 

15%-20% 

67.4 

± 

4.3 

144.8 

± 

23.8 

66.5 

± 

4.0 

1.97 

± 

0.17 

0.92 

± 

0.16 

20%-25% 

56.6 

± 

4.4 

118.7 

± 

11.5 

55.6 

± 

3.3 

1.96 

± 

0.19 

0.94 

± 

0.11 

25%-30% 

48.7 

± 

4.9 

100.0 

± 

12.0 

46.4 

± 

2.8 

1.91 

± 

0.22 

0.93 

± 

0.12 

30%-35% 

40.4 

± 

4.5 

81.1 : 

± 

10.4 

38.6 

± 

2.3 

1.91 

± 

0.24 

0.95 

± 

0.13 

35%-40% 

32.3 

± 

4.1 

63.3 

± 

6.1 

32.0 

± 

1.9 

1.98 

± 

0.28 

1.01 

± 

0.11 


TABLE XXVI. Charged particle multiplicity results for 62.4 GeV Cu+Cu collisions. The uncertainties include the total 
statistical and systematic uncertainties. 

Centrality_ (Aipart )_ {Ngp ) _ dNch/dr] _ (dVch/d??)/(0.5Vpart) _ (dNch/drj)/(0.5Nqp) 


0%-5% 

100.5 

± 

4.5 

229.3 ± 8.5 

135.3 

± 

11.1 

2.69 

± 

0.25 

1.18 

± 

0.11 

5%-10% 

88.3 

± 

4.8 

197.8 ± 15.0 

116.6 

± 

9.5 

2.64 

± 

0.26 

1.18 

± 

0.13 

10%-15% 

78.2 

± 

4.3 

171.7 ± 25.2 

97.8 

± 

8.0 

2.50 

± 

0.25 

1.14 

± 

0.19 

15%-20% 

67.4 

± 

4.3 

144.8 ± 23.8 

81.0 

± 

6.6 

2.40 

± 

0.25 

1.12 

± 

0.20 

20%-25% 

56.6 

± 

4.4 

118.7 ± 11.5 

67.0 

± 

5.5 

2.37 

± 

0.27 

1.13 

± 

0.14 

25%-30% 

48.7 

± 

4.9 

100.0 ± 12.0 

55.3 

± 

4.5 

2.27 

± 

0.30 

1.11 

± 

0.16 

30%-35% 

40.4 

± 

4.5 

81.1 ± 10.4 

45.3 

± 

3.7 

2.24 

± 

0.31 

1.12 

± 

0.17 

35%-40% 

32.3 

± 

4.1 

63.3 ± 6.11 

36.9 

± 

3.0 

2.28 

± 

0.35 

1.17 

± 

0.15 


TABLE XXVII. Transverse energy results for 200 GeV Cu+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 


Centrality 

(Apart ) 

(Ngp) 

dEr/dr] [GeV] 

(dAT/dr?)/(0.5Apart) [GeV] 

{dET/d'n)/{0.5Ngp) [GeV] 

0%-5% 

189.0 ± 5.2 

463.8 ± 17.6 

288.3 ± 17.3 

3.05 ± 0.20 

1.24 ± 0.09 

5%-10% 

164.2 ± 4.3 

400.3 ± 14.8 

249.8 ± 15.0 

3.04 ± 0.20 

1.25 ± 0.09 

10%-15% 

142.4 ± 3.7 

341.7 ± 12.7 

212.8 ± 12.8 

2.99 ± 0.20 

1.25 ± 0.09 

15%-20% 

122.6 ± 3.3 

288.9 ± 10.7 

179.4 ± 10.8 

2.93 ± 0.19 

1.24 ± 0.09 

20%-25% 

104.5 ± 3.5 

240.5 ± 11.0 

150.0 ± 9.0 

2.87 ± 0.20 

1.25 ± 0.09 

25%-30% 

88.5 ± 4.0 

199.0 ± 11.8 

124.5 ± 7.5 

2.81 ± 0.21 

1.25 ± 0.11 

30%-35% 

73.8 ± 3.6 

162.6 ± 9.8 

102.3 ± 6.1 

2.77 ± 0.21 

1.26 ± 0.11 

35%-40% 

60.9 ± 3.6 

131.0 ± 8.8 

83.3 ± 5.0 

2.74 ± 0.23 

1.27 ± 0.11 

40%-45% 

49.7 ± 3.2 

103.4 ± 8.8 

67.0 ± 4.0 

2.69 ± 0.24 

1.29 ± 0.13 

45%-50% 

39.9 ± 3.1 

80.6 ± 8.5 

53.1 ± 3.2 

2.66 ± 0.26 

1.32 ± 0.16 

50%-55% 

31.4 ± 3.2 

62.3 ± 8.1 

41.4 ± 2.5 

2.64 ± 0.31 

1.33 ± 0.19 

55%-60% 

24.3 ± 2.8 

47.1 ± 6.9 

31.9 ± 1.9 

2.63 ± 0.34 

1.36 ± 0.21 
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TABLE XXVIII. Charged particle multiplicity results for 200 GeV Cu+Au collisions. The uncertainties include the total 
statistical and systematic uncertainties. 


Gentrality 

(A^part ) 

) 

dNch/dri 

{dNch/dl]) / (0.5iVpart ) 

{dNch/dij)/{0.5Nqp) 

0%-5% 

189.0 ± 5.2 

463.8 ± 17.6 

333.5 ± 25.0 

3.53 ± 0.28 

1.44 ± 0.12 

5%-10% 

164.2 ± 4.3 

400.3 ± 14.8 

288.0 ± 21.4 

3.51 ± 0.28 

1.44 ± 0.12 

10%-15% 

142.4 ± 3.7 

341.7 ± 12.7 

244.5 ± 18.1 

3.43 ± 0.27 

1.43 ± 0.11 

15%-20% 

122.6 ± 3.3 

288.9 ± 10.7 

205.4 ± 15.1 

3.35 ± 0.26 

1.42 ± 0.11 

20%-25% 

104.5 ± 3.5 

240.5 ± 11.0 

171.2 ± 12.5 

3.28 ± 0.26 

1.42 ± 0.12 

25%-30% 

88.5 ± 4.0 

199.0 ± 11.8 

141.5 ± 10.2 

3.20 ± 0.27 

1.42 ± 0.13 

30%-35% 

73.8 ± 3.6 

162.6 ± 9.8 

115.9 ± 8.3 

3.14 ± 0.27 

1.43 ± 0.13 

35%-40% 

60.9 ± 3.6 

131.0 ± 8.8 

94.0 ± 6.7 

3.09 ± 0.29 

1.43 ± 0.14 

40%-45% 

49.7 ± 3.2 

103.4 ± 8.8 

75.2 ± 5.4 

3.03 ± 0.29 

1.45 ± 0.16 

45%-50% 

39.9 ± 3.1 

80.6 ± 8.5 

59.5 ± 4.2 

2.98 ± 0.31 

1.48 ± 0.19 

50%-55% 

31.4 ± 3.2 

62.3 ± 8.1 

46.3 ± 3.3 

2.95 ± 0.37 

1.48 ± 0.22 

55%-60% 

24.3 ± 2.8 

47.1 ± 6.9 

35.4 ± 2.5 

2.91 ± 0.39 

1.50 ± 0.24 


TABLE XXIX. Transverse energy results for 193 GeV U+U collisions. The uncertainties include the total statistical and 
systematic errors. 


Gentrality 

(Afpart ) 

dEr/dr^ [GeV] 

{dET/dr,)/{0.5Np^rt) [GeV] 

0%-5% 

418.8 ± 5.0 

783.0 ± 46.1 

3.74 ± 0.22 

5%-10% 

353.2 ± 6.0 

625.6 ± 36.9 

3.54 ± 0.22 

10%-15% 

296.7 ± 6.1 

504.0 ± 29.7 

3.40 ± 0.21 

15%-20% 

248.9 ± 6.8 

406.2 ± 23.9 

3.26 ± 0.21 

20%-25% 

207.6 ± 6.7 

325.9 ± 19.2 

3.14 ± 0.21 

25%-30% 

172.5 ± 6.5 

259.2 ± 15.3 

3.00 ± 0.21 

30%-35% 

141.6 ± 6.8 

203.7 ± 12.0 

2.88 ± 0.22 

35%-40% 

114.9 ± 6.9 

157.8 ± 9.3 

2.75 ± 0.23 

40%-45% 

91.8 ± 6.4 

119.9 ± 7.1 

2.61 ± 0.24 

45%-50% 

72.0 ± 6.2 

89.16 ± 5.3 

2.48 ± 0.26 


TABLE XXX. 
and systematic 

Charged particle multiplicity results for 193 GeV U+U collisions. ' 
errors. 

The uncertainties include the total statistical 

Gentrality 


(Vpart ) 

dNch/dri 

{dNch/dr]) / {0.5Nps,rt) 

0%-5% 


418.8 ± 5.0 

830.4 ± 67.8 

3.97 

± 

0.33 

5%-10% 


353.2 ± 6.0 

689.2 ± 55.5 

3.90 

± 

0.32 

10%-15% 


296.7 ± 6.1 

565.5 ± 44.9 

3.81 

± 

0.31 

15%-20% 


248.9 ± 6.8 

459.6 ± 36.1 

3.69 

± 

0.31 

20%-25% 


207.6 ± 6.7 

369.7 ± 28.7 

3.56 

± 

0.30 

25%-30% 


172.5 ± 6.5 

293.9 ± 22.6 

3.41 

± 

0.29 

30%-35% 


141.6 ± 6.8 

230.6 ± 17.5 

3.26 

± 

0.29 

35%-40% 


114.9 ± 6.9 

178.1 ± 13.4 

3.10 

± 

0.30 

40%-45% 


91.8 ± 6.4 

135.0 ± 10.1 

2.94 

± 

0.30 

45%-50% 


72.0 ± 6.2 

100.0 ± 7.4 

2.78 

± 

0.32 
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TABLE XXXI. Transverse energy results for 200 GeV d+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 


Centrality 

(Apart ) 

{Ngp 

) 

dEr/dr] [GeV] 

(d£;T/dj7)/(0.5Apart) [GeV] 

{dET/dr,)/{0.5Ngp) [GeV] 

0%-5% 

17.8 ± 1.2 

27.2 

± 

2.3 

20.3 ± 1.7 

2.29 

± 

0.24 

1.39 

± 

0.16 

5%-10% 

15.6 ± 1.0 

24.7 

± 

2.0 

17.4 ± 1.5 

2.24 

± 

0.23 

1.33 

± 

0.15 

10%-20% 

14.1 ± 0.9 

22.9 

± 

1.8 

15.4 ± 1.3 

2.18 

± 

0.23 

1.27 

± 

0.14 

20%-30% 

11.9 ± 0.7 

20.0 

± 

1.5 

13.2 ± 1.1 

2.21 

± 

0.22 

1.27 

± 

0.14 

30%-40% 

10.5 ± 0.6 

18.0 

± 

1.2 

11.3 ± 0.9 

2.16 

± 

0.22 

1.22 

± 

0.13 

40%-50% 

8.7 ± 0.5 

15.3 

± 

0.9 

9.5 ± 0.8 

2.20 

± 

0.22 

1.23 

± 

0.12 

50%-60% 

7.1 ± 0.4 

12.7 

± 

0.6 

7.8 ± 0.7 

2.19 

± 

0.23 

1.23 

± 

0.12 

60%-70% 

5.7 ± 0.4 

10.4 

± 

0.5 

6.3 ± 0.5 

2.21 

± 

0.23 

1.23 

± 

0.12 


TABLE XXXII. Charged particle multiplicity results for 200 GeV d+Au collisions. The uncertainties include the total statistical 
and systematic uncertainties. 


Centrality 

(Apart ) 

(A,p 

) 

dNch/dr; 

(dAch/d?7)/(0.5Apart) 

(dNch/drj)/ (O.SAqp) 

0%-5% 

17.8 ± 1.2 

27.2 

± 

2.3 

20.8 ± 1.5 

2.43 

± 

0.23 

1.53 

± 

0.17 

5%-10% 

15.6 ± 1.0 

24.7 

± 

2.0 

17.7 ± 1.2 

2.36 

± 

0.22 

1.43 

± 

0.15 

10%-20% 

14.1 ± 0.9 

22.9 

± 

1.8 

15.5 ± 1.1 

2.28 

± 

0.21 

1.35 

± 

0.14 

20%-30% 

11.9 ± 0.7 

20.0 

± 

1.5 

13.2 ± 0.9 

2.30 

± 

0.21 

1.32 

± 

0.13 

30%-40% 

10.5 ± 0.6 

18.0 

± 

1.2 

11.2 ± 0.8 

2.22 

± 

0.20 

1.25 

± 

0.12 

40%-50% 

8.7 ± 0.5 

15.3 

± 

0.9 

9.3 ± 0.7 

2.23 

± 

0.20 

1.22 

± 

0.11 

50%-60% 

7.1 ± 0.4 

12.7 

± 

0.6 

7.5 ± 0.5 

2.18 

± 

0.20 

1.18 

± 

0.10 

60%-70% 

5.7 ± 0.4 

10.4 

± 

0.5 

5.8 ± 0.4 

2.12 

± 

0.20 

1.12 

± 

0.10 


TABLE XXXIII. Transverse energy results for 200 GeV ®He+Au collisions. The uncertainties include the total statistical and 
systematic uncertainties. 

Centrality (Vpart )_ {N^p ) dEr/dri [GeV] {dET/drj)/{0.5Np^rt) [GeV] {dET/drj)/{0.5Ngp) [GeV] 


0%-5% 

25.0 

± 

2.0 

37.5 

± 

3.1 

26.7 

± 

1.8 

2.13 

± 

0.19 

1.33 

± 

0.14 

5%-10% 

22.6 

± 

1.3 

34.3 

± 

2.4 

23.2 

± 

1.5 

2.06 

± 

0.18 

1.27 

± 

0.12 

10%-20% 

19.9 

± 

1.1 

30.6 

± 

2.2 

20.6 

± 

1.4 

2.07 

± 

0.18 

1.26 

± 

0.12 

20%-30% 

17.0 

± 

1.0 

26.6 

± 

1.8 

17.7 

± 

1.2 

2.08 

± 

0.18 

1.25 

± 

0.11 

30%-40% 

13.8 

± 

0.7 

21.9 

± 

1.3 

14.9 

± 

1.0 

2.16 

± 

0.18 

1.27 

± 

0.11 

40%-50% 

10.9 

± 

0.7 

17.4 

± 

0.8 

12.0 

± 

0.8 

2.22 

± 

0.20 

1.29 

± 

0.10 

50%-60% 

8.16 

± 

0.5 

13.2 

± 

0.8 

9.3 

± 

0.6 

2.29 

± 

0.21 

1.32 

± 

0.12 

60%-70% 

6.01 

± 

0.4 

9.7 

± 

0.5 

7.0 

± 

0.5 

2.33 

± 

0.21 

1.34 

± 

0.11 


TABLE XXXIV. Charged particle multiplicity results for 200 GeV ®He+Au collisions. The uncertainties include the total 
statistical and systematic uncertainties. 

Centrality _(Vpart ) _ {Ngp ) _ dN^h/dr) _(dVch/d??)/(0.5jVpart)_ {dNch/dr])/ (O.SV^p) 


0%-5% 

25.0 

± 

2.0 

37.5 

± 

3.1 

26.3 ± 

1.8 

2.10 

± 

0.22 

1.40 

± 

0.15 

5%-10% 

22.6 

± 

1.3 

34.3 

± 

2.4 

22.7 ± 

1.6 

2.01 

± 

0.18 

1.32 

± 

0.13 

10%-20% 

19.9 

± 

1.1 

30.6 

± 

2.2 

19.9 ± 

1.4 

2.00 

± 

0.18 

1.30 

± 

0.13 

20%-30% 

17.0 

± 

1.0 

26.6 

± 

1.8 

16.9 ± 

1.2 

1.99 

± 

0.18 

1.27 

± 

0.12 

30%-40% 

13.8 

± 

0.8 

21.9 

± 

1.3 

14.0 ± 

1.0 

2.04 

± 

0.18 

1.28 

± 

0.12 

40%-50% 

10.9 

± 

0.7 

17.4 

± 

0.8 

11.2 ± 

0.8 

2.06 

± 

0.19 

1.28 

± 

0.11 

50%-60% 

8.16 

± 

0.5 

13.2 

± 

0.8 

8.4 ± 

0.6 

2.06 

± 

0.20 

1.27 

± 

0.12 

60%-70% 

6.01 

± 

0.4 

9.72 

± 

0.5 

5.9 ± 

0.4 

1.98 

± 

0.19 

1.22 

± 

0.11 
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